Finite size lattice results for the two-boundary 
Temperley-Lieb loop model 



Anita Kristine Ponsaing 

Submitted in total fulfilment of the requirements of 
the degree of Doctor of Philosophy 

June 2010 



Department of Mathematics and Statistics 
The University of Melbourne 



Abstract 



This thesis is concerned with aspects of the integrable Tempcrlcy-Lieb loop (TL(n)) 
model on a vertically infinite lattice with two non-trivial boundaries. The TL(n) model 
is central in the field of integrable lattice models, and diflrerent values of n relate to 
different physical models. For instance, the point n = relates to critical dense polymers 
and a corresponding logarithmic conformal field theory. The point n = 1 corresponds to 
critical bond percolation on the square lattice, and has connections with a combinatorial 
counting problem of alternating sign matrices and their symmetry classes. For general 
n, the TL(n) model is closely related to the XXZ quantum spin chain and the 6-vertex 
model. 

We construct the transfer matrix of the model, which describes the weights of all the 
possible configurations of one row of the lattice. When n = 1 the ground state eigenvec- 
tor of this matrix can be interpreted as a probability distribution of the possible states of 
the system. Because of special properties exhibited by the transfer matrix at n = 1, we 
can show that the eigenvector is a solution of the g-deformed Knizhnik-Zamolodchikov 
equation, and we use this fact to explicitly calculate some of the components of the eigen- 
vector. In addition, recursive properties of this transfer matrix allow us to compute the 
normalisation of the eigenvector, and show that it is the product of four Weyl characters 
of the symplectic group. Previous work in this area has produced results for the TL(1) 
loop model with periodic boundary conditions, two trivial boundaries and mixed (one 
trivial, one non-trivial) boundaries, but until recently little progress had been made on 
the case with two non-trivial boundaries. This boundary condition lends itself to calcu- 
lations relating to horizontal percolation, which is not possible with the other boundary 
conditions. 

One of these calculations is a type of correlation function that can be interpreted as 
the density of percolation cluster crossings between the two boundaries of the lattice. It 
is an example of a class of parafermionic observables recently introduced in an attempt 
to rigorously prove conformal invariance of the scaling limit of critical two-dimensional 
lattice models. It also corresponds to the spin current in the Chalker-Coddington model 
of the quantum spin Hall effect. We derive an exact expression for this correlation 
function, using properties of the transfer matrix of the TL(1) model, and find that it can 
be expressed in terms of the same symplectic characters as the normalisation. 

In order to better understand these solutions, we use Sklyanin's scheme to perform 
separation of variables on the symplectic character. We construct an invertible separating 
operator that transforms the multivariate character into a product of single variable 
polynomials. Analysing the asymptotics of these polynomials will lead, via the inverse 
transformation, to the asymptotic limit of the symplectic character, and thus to the 



asymptotic limit of the ground state normalisation and correlation function of the loop 
model. We construct the separating operator by viewing the symplectic characters as 
eigenfunctions of a quantum integrable system, and also explicitly construct the factorised 
Hamiltonian for this system. 
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Introduction 



Statistical mechanics is the study of systems composed of large numbers of small particles. 
The aim of the study is to investigate how global, collective effects arise from small- 
scale local interactions between the particles. The systems in question are too large to 
make accurate predictions on the level of individual particles, so the focus is instead on 
statistical quantities. Examples include the average density of particles and correlation 
functions between two chosen points. These statistical quantities can be used to calculate 
derived quantities such as heat capacity, energy densities and magnetisation. 

Global quantities of a system are usually expressed in terms of the system's temper- 
ature, as well as external parameters such as a magnetic field. As these parameters pass 
certain values, the global quantities can be observed to change dramatically, for instance 
at the boiling point of water. These values arc known as critical points. The location 
and nature of any critical points of a system are more examples of quantities of interest. 

Near a critical point, global quantities can be approximated by a power law (say, 
(T — Tc)", where Tc is the critical temperature). The exponent a in this approximation is 
called a critical exponent, and indicates the strength of the divergence as the parameter 
passes the critical point. It is widely believed (and experimentally verified) that systems 
with the same dimensionality, symmetries, and range of interactions will have the same 
critical exponents — a property known as universality. 

A well-known statistical mechanical model is for magnetisation in a bar of iron. Each 
individual atom can be thought of as a small magnet, whose alignment affects the align- 
ment of its neighbours. The strength of the effect depends on various factors, including 
the temperature of the bar and the external magnetic field. At low overall temperatures 
the atoms will prefer to keep their current alignment, but they will also try to align with 
their neighbours and the external field. This leads to the phenomenon of the bar having 
a spontaneous magnetisation if the external field is turned off. At high temperatures the 
particles are more likely to lose their alignment, so there is no spontaneous magnetisa- 
tion. The temperature at which the spontaneous magnetisation disappears is the critical 
temperature, and much interesting physics happens near this point. 

The local interactions of a system can be described by an energy function E. This 
can be variously interpreted as a function, a matrix, or more generally an operator. By 
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acting on a possible state s with the energy function, we get the energy E{s) associated 
with that state. When the system is in thermodynamic equiUbrium at a temperature T, 
the probability of observing a state decreases as the energy of that state increases. More 
precisely, the probability is 

1 E(s) 

P{s) = -e fcT , 

where k is Boltzmann's constant (approx. 1.38 x 10^"^^ J/ K), and the factor jg 
known as the Boltzmann weight of the state s. The normalisation Z is defined as 

s 

and is known as the partition function (for a derivation of these expressions see [87]). 
The partition function holds all the information about the probabilities of all possible 
states, and as such is an important quantity, especially for calculating expected values of 
observable properties. The partition function is also closely related to the thermodynamic 
free energy F of the system, by 

F = -ferinZ, Z = e-^. 

The free energy is the amount of energy in the system that is available to do work. 

The expected value of an observable A is the weighted sum over the values of A for 
all possible states, weighted by the probability of each state. In short, 

1 E(s) 

{A) = -Y,A{s)e-^ . 

s 

Another important physical quantity is the correlation function, which describes how two 
spatially separated points affect each other's behaviour. For a quantity A, the correlation 
function is given by 

(Ai^2)-(^l>(^2), 

where Ai is the quantity evaluated at the first point, and A2 is evaluated at the second. 

A system is said to be solvable if it is possible to analytically calculate the physical 
quantities described above. In other words, the partition function (or equivalently, the 
free energy) must be calculable [4]. For physical systems, this is usually an impossible 
task, and as such, most of the work done in this field is either based on numerical ap- 
proximations or simplified mathematical models. Nevertheless, these models can be quite 
powerful, as universality predicts that a model with the same dimensionality, symmetry 
and range of interactions as a physical system will also have exactly the same the critical 
exponents. 
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Lattice models 



Some of the most commonly used models in statistical mechanics are lattice models, 
where the particles are located at regular intervals, and interactions are defined between 
them. Calculating the partition function thus involves solving the combinatorial problem 
of finding all possible lattice configurations, and then summing the probabilities for each 
configuration; the latter is generally the more difficult task. 

Even though these models do not reflect reality in full detail, by letting the number of 
lattice sites tend to infinity they can be used to obtain a good idea of the behaviour of a 
physical system. This infinite limit is known as the thermodynamic limit, or equivalently 
the asymptotic limit. In addition, as previously mentioned, universality predicts that an 
appropriately designed lattice model will have the same critical exponents as a physical 
system. 

One of the most useful tools in the study of lattice models is the method of transfer 
matrices. Roughly speaking, a transfer matrix describes all possible configurations of one 
row of the lattice. If the height and width of the lattice are M and L respectively, it can 
also be shown that in the large M limit, the partition function goes like 

Zlm ~ Af^, 

where A^, is the largest eigenvalue of the transfer matrix (corresponding to the state with 
the lowest energy, or the ground state). The thermodynamic limit then corresponds to 
taking the limit as L goes to infinity. 

The transfer matrix can be constructed from operators called i?-matrices that act 
locally. If it can be shown that the i?-matrices satisfy the Yang-Baxter equation, then it 
is possible to construct an infinite family of commuting transfer matrices. This infinite 
family is parameterised by the so-called spectral parameter. The fact that the family 
of transfer matrices commute impUcs that they arc simultaneously diagonalisabic and 
share eigenvectors, and since each member of the family has a dift'erent value of the 
spectral parameter, these eigenvectors cannot depend on the spectral parameter. This 
fact provides access to a host of methods for finding the eigenvalues and eigenvectors of 
the transfer matrix, and often leads to the system being solvable. 

The most famous of these methods is the Bethe ansatz method. This method works 
by guessing a form for the eigenvectors of the transfer matrix (called the ansatz). The 
ansatz depends on a set of parameters, which are fixed by constraint equations resulting 
from the transfer matrix eigenvalue equation. The Bethe ansatz method relies on the 
existence of a pseudovacuum — an eigenvector of the transfer matrix for which there 
exists an annihilation operator. 

Perhaps the most widely known lattice model is the Ising model, which is a model of a 
magnet, as described earlier. The Ising model was posed by Lenz in 1920, and solved for 
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one dimension by Ising in 1925 [45]. In two dimensions the solution for zero external field 
was given by Onsager in 1944 [74]. For non-zero external magnetic field, Zamolodchikov 
found a solution in 1989 for a field-theoretic version of the two-dimensional Ising model in 
the thermodynamic limit [108]. In 1992, Warnaar, Nienhuis and Seaton found a solution 
for a two-dimensional model in the same universality class [102]. 

The Ising model is typically defined on a square lattice. At each site is a molecule 
with a magnetic moment (Tj, which can take the value -|-1 (up) or —1 (down). These 
magnetic moments can also be regarded as spins. Each spin interacts only with its 
closest neighbours, and the interaction depends on the values of the spins and a coupling 
constant J . The energy function of the Ising model is 

<ij> i 

where the first sum is over all nearest neighbours, and B is the strength of the external 
field. When JT" < 0, this system models a ferromagnet, since the energy is lower when 
adjacent spins are aligned. When > the system models an antiferromagnet. In two 
dimensions, the model exhibits spontaneous magnetisation below a critical temperature, 
but in one dimension no such behaviour occurs. 

Naturally this is a highly idealised model, as real magnetic moments can point in 
more than two directions. Several generalisations of the Ising model have since been 
proposed, including the Potts model and the n-vector model. 

The Q-state Potts model, or simply Potts model [4,101,104], has scalar spins that 
can take one of Q values, and the interactions between nearest neighbours are given by 
Kronecker delta functions. This is a well-known and widely studied model in the field of 
two-dimensional lattice models. In its current form it was first written down by Potts in 
1952 [79]. The energy function of the Potts model is 

<ij> 

and the partition function is then 

a 

where the sum is over all possible spin configurations. This can be rewritten as 

o- <ij> 

with u = ekT —1. 
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The Potts model can be reformulated as a loop model in the following way [4]. Con- 
sider the system of spins as a planar graph L where the spins sit on vertices and edges 
are drawn between neighbouring spins. On each edge, a bond may be drawn according to 
the following rules: if a pair of neighbouring vertices have different spin values, no bond 
is drawn; if the vertices share the same spin value, there is a choice of drawing a bond 
on the edge between them (with a weight of n), or leaving the edge blank (with weight 
1). We refer to an arrangement of bonds on £ as a bond-graph. Then the expansion of 
each product in the partition function represents all the possible bond-graphs associated 
with a spin configuration. An example of this is shown in Figure 1. 



Figure 1: Possible bond-graphs for the spin configuration ((1,2), (1,1)). 

By considering instead the spin configurations possible on a particular bond-graph, 
we can turn the partition function into a sum over bond-graphs. Let C be the number 
of connected components of the bond-graph (clusters of bonds, including single vertices). 
All the vertices in a connected component must share a spin value. This spin value could 
be any of the Q, and since each connected component is independent of the others, there 
are ways to choose these values for the entire bond-graph. A bond-graph with I bonds 
has a weight of u'. Then the partition function is a sum over all the possible bond-graphs 
G on 

G 

The connection to the loop model is made by drawing loops around the bonds in the 
following two ways: 

...).(... ^ 

Every cluster has Q possible spin values, which leads to every closed loop in the loop 
expansion having a weight of \fQ. 

At Q = 1 the loop model is equivalent to the two-dimensional square lattice bond 
percolation model [77, 90] . In this model, bonds are drawn between neighbouring sites 
of the lattice with probability ■p. On an infinite lattice, percolation occurs when there is 
an infinite sized cluster of bonds. This model has a second order phase transition (that 
is, the first derivative of the free energy is continuous but the second is not) at critical 
probability p = 1/2, above which percolation occurs with probability 1. 

By orienting the loops, the Potts loop model can be generalised to a directed loop 
model, which can itself be mapped to the six vertex model. These mappings are explained 
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in more detail in Chapter 12 of [4], and in [5]. 

The six vertex model is another classical two-dimensional lattice model, also known 
as the square ice model. On a square lattice, oxygen atoms are placed at each vertex, 
and on the edges between them are placed hydrogen atoms. A single water molecule 
consists of one oxygen bonded with two hydrogens, and correspondingly the 'ice rule' 
for the lattice model states that of the four hydrogens surrounding an oxygen, two will 
be near the oxygen and two will be distant from it. The ice rule leads to six possible 
configurations around a vertex: 

In H -5 „ 

H 

H OH HO H H O H 

H H H 



H 

H-6 H 



H 

H 6-H 



6 H 
H-O-H 
H 



The six possibilities arc usually drawn as arrow configurations, where each arrow cor- 
responds to a hydrogen atom and points to the nearest oxygen. Each vertex configuration 




Figure 2: The 6 vertex configurations with associated Boltzmann weights a, b, 
and c. 



has a Boltzmann weight associated with it as shown in Figure 2. The Boltzmann weights 
a, b, and c appear in the partition function of the model, 

conf. 

where the sum is over all possible configurations of the lattice and nia, rrib and nic are 
the numbers of vertices of types 1 &: 2, 3 &; 4, and 5 & 6 respectively. 

The six vertex model has strong connections to combinatorics. When defined on a 
finite MxM grid, with 'domain wall' boundary conditions — outward-pointing arrows at 
the top and bottom, and inward-pointing arrows at the sides — the possible configurations 
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of the lattice are in bijection with MxM alternating sign matrices (ASMs). The bijection 
is made by replacing all vertices of type 1, 2, 3, and 4 with a 0, vertices of type 5 with a 
— 1 and type 6 with a 1. 



1 

1 -1 1 

1 



A similar bijection can also be made with fully packed loop diagrams (FPLs). For a 
thorough review see [80] . 

The six vertex model was solved (that is, the partition function was calculated) in 
1967 by Lieb [60,61] and Sutherland [98]. Lieb also found that the transfer matrix of 
the six vertex model shares eigenvectors with the Hamiltonian of the one-dimensional 
quantum XXZ spin chain, which was solved by Yang and Yang using the Bethe ansatz 
in the previous year [105, 106]. 

The XXZ model is a quantum generalisation of the Ising model. Its Hamiltonian is a 
sum over nearest neighbours, 



where a^, cr^, and are spin operators represented by the Pauli matrices, and A is 
known as the anisotropy parameter. When A = 1, the system becomes the Heisenberg 
model, and the value A = corresponds to the free fermion point, where the model is 
related to the Ising model. 

The eigenvectors of the six vertex transfer matrix only depend on the Boltzmann 
weights in a certain combination, which, as it turns out, exactly corresponds to the XXZ 
anisotropy parameter. 



A = 



+ 6^ - 
2ab 



In fact, the transfer matrix of the six vertex model and the Hamiltonian of the one- 
dimensional quantum XXZ spin chain are related through the Taylor expansion of the 
transfer matrix, 

TQy{n) = T6v(0) + n Hxxz + ...), 

where /x is the spectral parameter. 

This relationship is an example of a common occurrence in classical two-dimensional 
lattice models, where the logarithmic derivative of the transfer matrix taken at /x = 
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corresponds to the Hamiltonian of a related one-dimensional quantum system. The 
Hamiltonian is usually a local operator, unlike the transfer matrix, so it is often easier 
to deal with, however it can also contain less information. 

Another generalisation of the Ising model is the n-vector model (also known as the 
0(n) model), introduced by Stanley in 1968 [97]. In this model, each spin is a unit vector 
with n components. The energy function of the model is 

<ij> 

The name 0(n) refers to the fact that the spins have the symmetry of the orthogonal 
group. The limit n — > is related to the theory of self-avoiding walks, which has an 
application in the study of polymers. 

The n-vector model, similarly to the Potts model, has a reformulation as a loop model 
(see [7,8,35,55,68]), and this allows for a generalisation such that n need not be a positive 
integer. The 0(n) loop model on a square lattice has nine possible states for each lattice 
face: 




The loops decorating the faces (or plaquettes) correspond to interactions between spins 
(the spins are located on the centres of the edges) . Each plaquette has a weight associated 
with it, and various results have been obtained for systems with different restrictions on 
the weights. 

The completely packed version [66] , where the weights of all but the last two pictures 
are set to zero, is the one this thesis will be focusing on. The model has a rich underlying 
mathematical structure, best understood in terms of the Temperley-Lieb algebra. For 
this reason we will refer to the completely packed 0(n) loop model as the Temperley-Lieb 
loop model, or the TL(n) model. This model is related to the Q-state Potts model on 
the square lattice by ^/Q = n. 

The Temperley-Lieb loop model 



It would seem natural to define the Tcmpcrlcy-Lieb loop model in terms of its origins 
as the n-vector model. However the connection is mathematically involved [8, 33, 35, 55] 
and not altogether intuitive. Further, in recent years the TL(n) model has garnered a 
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significant amount of interest quite separate from the spin model formulation. For these 
reasons we will define it directly as a loop model. 

The TL(n) loop model is defined on a horizontally finite, vertically semi-infinite square 
lattice, where each face of the lattice has loops drawn on it in one of two orientations, as 
described earlier. We choose the infinite direction to be downwards, but this is simply a 
convention. Weights are then given to closed loops, and we are interested in the remaining 
connectivities at the top of the lattice. The lattice may be wrapped around a cylinder 
to form the periodic TL(n) model, or placed on a finite width strip with various types of 
boundary conditions at the left and right. If loops are allowed to end on a boundary, we 
call the boundary open (or non-diagonal), otherwise the boundary is closed (or diagonal). 
In addition to the periodic model, there arc three notable versions of the model; namely 
the zero-boundary model with closed boundaries at the left and right, the one-boundary 
model with one closed and one open boundary, and the two-boundary model with two 
open boundaries. We are interested in the latter case. 

The transfer matrix of the TL(n) model describes all the possible configurations on 
one row of the semi-infinite lattice, and one can add rows to the lattice by acting with the 
transfer matrix at the top of the semi-infinite strip. When n = 1 the largest eigenvalue of 
the transfer matrix is 1, so the partition function is trivial.* These properties allow us to 
calculate other statistical quantities of interest for the two-boundary model with n = 1, 
such as the probability of observing a certain configuration at the top of the lattice, and 
a correlation function between two spatially separated points. 

The TL(n) model, at special values of n, is closely related to a number of different 
statistical mechanical models. At n = 1 there is an easy mapping to the two-dimensional 
bond percolation problem, which is made by drawing hulls around the percolation clusters 
as pictured: 



Percolation of the bonds then corresponds to percolation of the hulls, or loops. Also 
at n = 1, the action of the transfer matrix on the lattice of the TL(1) model can be 
interpreted as the evolution of a stochastic process (the stochastic Raise and Peel model 
[78]), and the ground state eigenvector of the transfer matrix is the steady state of 
this process. In addition, the TL(ri) model is closely related to the quantum XXZ spin 

*This fact is in agreement with the triviahty of the Potts model at Q = 1, when every spin has the 
same value. There exists a sypersymmetry interpretation of this point, see [40]. 
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chain in one dimension, as the Hamiltonians of the two models correspond to different 
representations of the underlying Temperley-Lieb algebra. The anisotropy parameter A 
in the XXZ model is related to n by n = — 2A. 

The ground state of the TL(1) loop model with periodic boundary conditions has also 
been connected with the combinatorial counting problems of alternating sign matrices, 
fully packed loops and plane partitions. This connection has been known as the Razumov- 
Stroganov (RS) conjecture [30,82-85], but was recently proved using combinatorial means 
by Cantini and Sportiello in [11]. The original conjecture also sparked a host of related 
conjectures, largely concerning different boundary conditions or topologies on the TL(1) 
model and different boundary conditions or symmetry requirements on the counting 
problems, which at the time of writing remain unproved [3,17,67]. Notable by its absence 
is a conjecture relating to the ground state of the TL(1) loop model with two distinct 
boundaries, which as yet has no connection to any similar combinatorial problem, though 
if the boundaries are identified a connection can be made [24]. 

A key element in the solutions of TL(n) model with different boundary conditions is 
the g-deformed Knizhnik-Zamolodchikov (gKZ) equation [10,23,26-28,30,48,49,75,109]. 
With the relationship n = —{q + q~^), it can be shown that the ground state of the TL(1) 
model satisfies the gKZ equation for q = e^^^^^. The version of the gKZ equation used 
depends on the boundary conditions. This connection provides a tool set for solving for 
the ground state, which would otherwise be inaccessible. 

Recently, a further connection has been made between the TL(1) loop model and 
the quantum spin Hall effect [12,43,52]. In particular the 2-boundary version is closely 
connected to this effect, since a type of correlation function between the two boundaries 
(see Chapter 3) corresponds to the spin current in a model of the quantum spin Hall 
effect. 

Layout of the thesis 

This thesis is concerned with properties of the TL(n) loop model with two open bound- 
aries, at the special point n = 1. Each horizontal position i on the lattice has a variable 
Zi associated with it, and the left and right boundaries have associated variables and 
respectively. In addition, the transfer matrix T depends on the spectral parameter w, 
and since it can be shown that [T(w),T{u)] = 0, the eigenvectors do not depend on this 
parameter. 

Chapter 1 defines the TL(n) loop model with general parameters, first by defining 
the Temperley-Lieb algebra and then describing the lattice model. The transfer matrix 
and Hamiltonian are defined, and the specialisation to n = 1 is described. 

The main results of Chapter 2 are the exact calculation of key components of the 
ground state eigenvector of the 2-boundary transfer matrix for finite system size L, as 
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well as the normalisation of the eigenvector. These quantities are expressed in terms of 
the polynomial character of the symplectic group, which is the Schur polynomial for the 
root system of type C. The reason why the remaining components of the eigenvector 
become harder and harder to specify as the system size increases is discussed. 

Another quantity that is exactly calculated for finite system size, in Chapter 3, is the 
boundary-to-boundary correlation function, which is equivalent to the average density of 
percolating clusters between two points in a two-dimensional percolation model, defined 
on an infinite lattice of finite width. This quantity can also be interpreted as the spin cur- 
rent in the Chalker-Coddington model of the quantum spin Hall effect. The correlation 
function is expressed in terms of the same symplectic characters as the normalisation. 

A key aspect of any statistical mechanical model is the thermodynamic limit. The 
results in Chapter 2 and Chapter 3 are obtained exactly for finite system sizes, which is 
usually not possible for critical systems. Generally it is only possible to derive approxi- 
mate results in the thermodynamic limit. Our calculations thus allow for the calculation 
of exact thermodynamic quantities, by taking the limit as L ^ oo. It will then be possible 
to compare our results to the previously calculated non-rigorous continuum limit [12]. 

In order to obtain the thermodynamic limit, we need to consider the large L asymp- 
totics of the symplectic character that appears in all the key results. In particular we are 
interested in the effect the two boundaries have in the asymptotic limit. It is therefore 
desirable to analyse the asymptotics of the symplectic character when all its arguments 
are set to 1, except for the two variables that correspond to the boundaries. In Chapter 4, 
a separation of variables method initiated by Sklyanin [92] is extended to the symplectic 
character. It is hoped that an asymptotic approximation for the separated polynomial 
will be easier to find than the symplectic character itself, and by reversing the process 
of separation the result can be used to calculate the thermodynamic limit of the phys- 
ical quantities associated with the TL(n) model. At present, however, this is an open 
problem. 
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Chapter 1 

The Temperley— Lieb loop model 



Sklyanin's transfer matrix [92] of the TL(n) model, and consequently the Hamiltonian, 
can be expressed in terms of algebraic generators satisfying a Temperley-Lieb algebra 
[100], see for example [27,77]. The imposed boundary conditions decide which particular 
version of the Temperley-Lieb (TL) algebra is needed. In this chapter we will introduce 
the TL(n) model on a strip with open boundaries on both sides, which can be described 
in terms of the two-boundary Temperley-Lieb (2BTL) algebra [18]. Models with two 
reflecting or diagonal boundaries, as well as with mixed boundaries were studied in [23, 
27,109]. The periodic version was considered in [66,82,84]. 

1.1 The Temperley-Lieb algebra 

The Temperley-Lieb algebra is an 5^2 quotient of the Hecke algebra (see Appendix A). 
The different types of the Hecke algebra lead to different boundary conditions of the 
TL algebra. The 2BTL algebra comes from the Hecke algebra of type BC [18,37], see 
also [16]. 

One of two distinguished representations of the 2BTL algebra [18] is in a space of 
connectivities or link patterns, described in Section 1.1.2. This representation is relevant 
for the TL(n) model with open boundaries [46,67], and we will use it to illustrate the 
generators and relations that define the TL algebras. 

We will give brief definitions of the TL algebra with the four most common boundary 
conditions, before describing the two-boundary version in more detail. We then use the 
2BTL algebra to define the TL(n) loop model for general n, and in Section 1.3 discuss 
the specialisation n = 1 at which we obtain physical results for the system. 
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Trivial boundaries 



Definition 1.1.1. The best known version of the TL algebra is generated by elements 
Cj, l<i<i^— 1, which satisfy relations 



e^ = —2 cos 7 Cj, 



Vi, 

1 < i < L - 2, 

2 < i < L - 1, 



where L is the system size and 7 is a complex parameter, related to the parameter n by 
n = —2 cos 7. 

In the link pattern representation, the generators look like strings across a strip, 

1 ... i i+l ... L 

ei = 



and multiplication in the algebra from the left corresponds to vertical concatenation of 
the pictures from the top. The relations can be put in terms of two rules: firstly, from 
the rule = —2 cos 7 Cj, closed loops are removed and replaced with a factor of —2 cos 7, 



i i+l 

o 



-2 cos 7 



! i+l 



and secondly, from the rule eiei±iei = Cj, strings are pulled tight, 
1 ... i-1 i i+l ... L 



TIT" 



1 



i i+l 
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Periodic boundaries 

Definition 1.1.2. The periodic version has the generators and relations from Defini- 
tion 1.1.1, as well as an extra generator that satisfies the relations [50,62] 

e\ = -2 cos 7 CL, 

\ k = l,L-l, 
eiCi = CjCL, 2<i<L-l. 

This algebra is infinite-dimensional, so when L is even an additional relation can be 
imposed, 

(eiea . . . cl-i) (6264 • • • e^) (eica . . . cl-i) = v (eics . . . cl-i). (l-l-l) 

This relation produces a finite-dimensional quotient of the periodic TL algebra. 

The pictorial representation is similar to the trivial boundary version, but this time 
the strip closes on itself to form the surface of a cylinder. 




In this representation, the relation ejCi-iei = Cj looks like 
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The quotient relation (1.1.1) describes the weight given to non-contractible loops, 





One boundary 

Definition 1.1.3. The one-boundary Tempcrley-Lieb (IBTL) algebra, otherwise known 
as the blob algebra [62], also has the generators and relations of the trivial boundary 
version Definition 1.1.1, along with a left boundary generator eo, depicted as 



eo 



sm cjo 



which satisfies 

eieoei = eo, = r^eo, 

sin(a;o + 7) 

where ojq is an additional complex parameter. Note that we do not impose eoeico = eg. 
Pictorially the first relation is 



1 2 



1 2 



and the second is 



) 



— sin ojQ 



sin(a;o + 7) K 



Both of these relations involve a closed loop at the left boundary. However, only the 
second relation produces a non-trivial factor. To make clear the distinction, we introduce 



15 



the notion of parity by supposing the generators are shaded in the following way: 
Ci = i even, Cj = , i odd, 

and eo has the same shading as the even case. Then a loop that has both ends connected 
to the left boundary produces a non-trivial factor iff the inside of the loop is shaded. 

Two boundeiries 

Finally, the version we will be working with is the 2BTL algebra [34], which has the 
generators of the one-boundary case above, as well as a right boundary generator sl 
depicted as 



which satisfies similar relations to eg, 



n 



2 —smujL 

CL-ieieL-i = CL, = -— ■ — tCl. 

sm{ujL + 7) 

Definition 1.1.4. Written all together, the algebraic relations for the 2BTL are 

ef = —2 cos 7 Ci, 1 < i < L — 1, 

e? = ^^^^e,, i = 0,L, (1.1.2) 

sm[uii + 7) 

eiCiiiej = Cj, 1 < i < L - 1, 

along with the idempotent relations (1.1.4) described in the next section. 

Loops connected to the left boundary behave as in the IBTL, but at the right bound- 
ary the parity of L must also be taken into account. If L is even, a loop that has both 
ends connected to the right boundary produces a non-trivial factor iff the inside of the 
loop is shaded, but if L is odd the factor appears iff the inside of the loop is not shaded. 



1.1.1 The idempotent relations 

The versions of the TL algebra with and 1 boundaries are both finite-dimensional. 
However, the periodic and 2-boundary versions are not, since it is possible to form a 
string winding all the way around the cylinder (for the periodic case) or stretching from 
the left boundary to the right (for the two-boundary case), as in Figure 1.1. 

It was shown in [18] that all finite-dimensional irreducible representations of the 2BTL 
satisfy two additional relations, which we will describe now. 
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Figure 1.1: The arrangement 62656163606264 when L = 5 has a loop stretching 
from the left boundary to the right boundary. 



Definition 1.1.5. We first define two (unnormalised) idempotents /i and I2 as follows: 



6163 •••ei-i, L even, 
6163- --eL, L odd, 



For example, when L = 6 the idempotents are 
1 2 3 4 5 6 

h= ^ 



6062 •••ei, L even, 

6062 •••eL_i, L odd. 



1 2 3 4 5 6 



(1.1.3) 



r\ 



The double quotient of the 2BTL algebra has the additional relations: 



hhh = bh, 



(1.1.4) 



where b is an additional parameter. In our pictorial representation, these relations have 
the effect of removing pairs of loops stretching from the left boundary to the right and 
replacing them with a factor of b. 

1.1.2 Link pattern space 

Here we will define the space LP^, spanned by link patterns (sometimes called connec- 
tivities) in terms of the loop representation of the 2BTL generators above. This space is 
equivalent to the space spanned by (a variant of) anchored cross paths [81]. The space 
LFl forms the state space of the TL(ra) loop model [91]. An example of a link pattern is 
given in Figure 1.2. 

Definition 1.1.6. A link pattern is a non-crossing matching of the integers 0, 1, ... , L+1. 

The matching between the integers 1, . . . ,L is pairwise, whereas and L + 1 may be 
matched with, or connected to, an arbitrary number of other integers. The integers 
and L + 1 are respectively referred to as the left and right boundary. 

We can express elements of LP^ as words in the generators of the 2BTL. We choose 
one of the idempotents to be the shortest word, and act with any combination of the 
2BTL generators. Then by the relations (1.1.2) and (1.1.4), the resulting picture will 
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12 ... 7 8 




Figure 1.2: A link pattern for L = 8. 



reduce to one of the link patterns of size L, multiplied by the weight as introduced by 
the relations. For example, with L = 2 we define the shortest word to be I\ = ei, and 
acting with the combination eoei produces 

eoci/i = —2 cos 7 cqIi. 



In fact, the four basis elements of LP2 are Ii, eo/i, 62/1 and 6062/1, and these are 
represented respectively by the pictures 

12 12 12 12 



Note the line connecting the two boundaries in the final picture. If we had defined 
the shortest word to be I2 = 6062 instead of /i, a line would appear on the first three 
pictures and not the final one. Using the word representation of link patterns and starting 
from Ii, the link pattern with 1 connected to the left boundary and 2 to the right will 
always have one of these lines, and the others never will. Because of this, and because it 
is always possible to find out whether a picture should have a line, we omit single lines 
connecting the two boundaries when referring to a link pattern. In this way, the above 
pictures can be represented as the link patterns in Figure 1.3. 



12 12 12 12 




Figure 1.3: The link patterns for L = 2, corresponding to h, eoh, 62/1 and 
6062/1 respectively. 

We use a shorthand notation for the link patterns, in terms of a sequence of opening 
'(' and closing ')' parentheses. The ith bracket in the sequence refers to whether the 
ith site is connected to some place to the right of it (opening bracket) or to the left 
(closing bracket). Since the loops arc non-crossing and every site from 1 to L must be 
connected, this provides a unique labelling. The above link patterns for L = 2 are thus 
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given respectively by 

0, )), ((, )(• 

Because we can independently place an opening or closing parenthesis at each site, the 
dimension of the space LP l of link patterns of size L for the two-boundary Temperley- 
Lieb algebra is 

dimLPi = 2^. 

1.1.3 Path representation 

We now present a representation of the link pattern space using the graphical depiction of 
Cj as a tilted square tile decorated with small loop segments; and of eg and as similarly 
decorated half-tiles lying against the boundaries of the picture, as shown in Figure 1.4. 



7 i i+1 L 




(a) eo (b) Ci (c) cl 



Figure 1.4: Tiles for the 2BTL generators. 

The shortest word, given by one of the two idempotents, is depicted as a row of tilted 
half-tiles, which lie along the bottom of the picture, as in Figure 1.5. 




1 2 3 4 5 6 7 



Figure 1.5: The idempotent Ji for L = 7. The label i on the x-axis refers to 
the centre of the tile e,. 

Multiplication in the 2BTL algebra corresponds to vertical concatenation of the tiles. 
As an example, below are shown the algebraic relations e| = —2 cos 7 Cj and SiCi-iei = e^. 
The other relations are similar. 

I i+l i i+1 i-1 i i+1 i-1 / i+1 

= —2 cos 7 
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Because of these relations, each Unk pattern corresponds uniquely to a path, which is 
traced out by the top of the tiles. A path is defined in this representation as beginning 
at the left boundary, ending at the right, and taking one step either up or down for each 
step to the right. The shorthand notation defined in the previous section works with 
the path representation as well, where steps down and up are symbolised by ')' and '(' 
respectively. It is easily seen (as in Figure 1.6) that when the path steps up, the loop 
in the link pattern connects somewhere to the right, and the opposite is true for down 
steps. 




1.1.4 Spin chain representation 

The 2BTL algebra has another distinguished representation, also of dimension 2^ [67]. 
This representation is in the tensor product space (C^)®^, giving rise to the quantum 
XXZ spin chain with non-diagonal boundary conditions on both sides. We define the 
Pauli spin matrices. 



a 



1 

1 

-i 

1 

1 

-1 
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and will use the notation 



^x,y,z ^ ^^i-l ^ ^x,y,z ^ ^®L-i_ 



In the spin chain representation the Temperley-Lieb generators take the form [19] 

ei^l {(Tfc^f+i + ^t^f+i + cos 7 - l) - i sin 7 {at - C7|+i)) , 

Co 1-^ r (cos^i cr? + sin^i + icoswo cr? — sinwo) , 

2 sm(wo + 7) 



(cos ^2 ctl + sin ^2 — icosw^ — sinw^) , 



2sin(a;L + 7) 

where 9i and 62 are related to the idempotent parameter b by 

cos(6'i - 62) + cos(7 + 1:^0 + ool) 



b= < 



2 sin(7 + Wo) sin(7 + lol) 

cos(gi - 02) + cos(a;o - lol) 
2sin(7 + a;o)sin(7 + a;L) ' 



L even, 
Lodd. 



Note that b only depends on the difference of these two parameters. Due to the rotational 
symmetry in the spin x-y plane, the extra freedom (^1 + 62) can be interpreted as a free 
gauge parameter on which the algebra does not depend. 

1.2 The TL(n) lattice 

The two-boundary TL(n) loop model is defined on a lattice of width L, vertically infinite 
with a boundary on the left and right hand sides. Each face of the lattice is decorated 
with loops, which can end on the boundary or close back on themselves. Each bulk 
(square) face has two possible configurations of loops. 



and 



The boundary (triangle) faces also have two possible configurations, one joining the loops 
on two adjacent rows of the lattice, the other connecting them both to the boundary. At 
the left boundary the possibilities are 



and 



and at the right boundary the above pictures are reflected. Because of the action of the 
boundary faces, it is easy to see that the rows of the lattice come in pairs, as illustrated 
in Figure 1.7. 
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Figure 1.7: One of the configurations of four rows of the lattice for L = 7. 

One of the important aspects of this model is how the loops connect different parts 
of the lattice. If we choose a configuration for each of the faces, and cut the lattice 
along a horizontal line between two pairs of rows, it is clear that below the line there 
arc a collection of open loop segments, which connect different horizontal positions (see 
Figure 1.8). In the example, position 1 is connected to the left boundary, position 2 to 5, 



Figure 1.8: The link pattern below the line is produced by the bottom two rows 
of the lattice in Figure 1.7. 

3 to 4, and 6 to 7. This is a clear example of a link pattern as defined in Section 1.1.2. As 
in the 2BTL algebra, we can ignore the specific paths the loops take and replace closed 
loops and loops that have both ends connected to a boundary by their respective weights. 
Once this is done the only thing left to consider is the remaining link pattern, which we 
call a. Every configuration of the semi-infinite lattice below the horizontal line has one 
of the 2^ link patterns at the top. 

The two rows of the lattice above the link pattern can have one of 2^-'"+^ possible 
configurations, and each one will map a link pattern into another. As an example, the 
configuration in Figure 1.8 produces the link pattern 

I ^ ^ I 
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and introduces the weight 4cos^7 from the two closed loops.* We represent the two 
lattice rows as a matrix T (called the transfer matrix) , with the ijih entry being the sum 
of weights produced by all the configurations that map the jth link pattern to the ith 
(where there is some ordering on the link patterns). 

A vector with a basis in the space of link patterns can be written as 

If this vector is the unique 'ground state' eigenvector of T, 

T|*) =A„,ax|*), 

where Amax is the maximum eigenvalue, then repeated action of T on some initial state 
I in) produces 

T*^|in) ~ A^^^l^r), as A: ^ oo. 

This process can be seen as building up a semi-infinite lattice with infinite copies of the 
transfer matrix. The eigenvector then expresses the relative weights of all possible link 
patterns of size L, produced by the configurations on the lattice. 

As shown by Di Francesco and Zinn-Justin for other types of boundary conditions 
[27,109], it is possible to derive exact closed form expressions for certain properties of the 
ground state eigenvector for finite system sizes. To achieve this one needs to generalise 
the model just described, by considering an inhomogeneous version of the transfer matrix. 
This is done in the next section. 

1.2.1 Baxterisation 

In order to define the transfer matrix wc will first introduce the operators R and K, as 
well as their unchecked versions, to be defined shortly. We furthermore list some useful 
properties that we will need in later calculations. Throughout the following we will use 
the notation [z] for 

[z\ = z- z"^ 

and define 

g = e^^. 

Definition 1.2.1. The Baxterised elements Ri{z), and the boundary Baxterised elements 

*The loop at the left boundary does not contribute a weight, as it only involves one left boundary tile, 
and is therefore a loop produced by the rule eieoei = ei. At the end of this chapter all the weights will 
be set to 1 and this technicality will no longer be important. 
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Kq{zX) and Kl{z,C,) of the Temperley-Lieb algebras are defined as 

1 [q/z] - [z] ei 



Ri{z) 



K,{z) [qz] 



J^i{z,Q)= ri—j — T 7 , i = 0,L, 

k{l/zX,i^i) 



(1.2.1) 



where the parameter z is called the spectral parameter. Each boundary element can also 
be equipped with an additional free parameter ^. 



The weight k{z, C,0Ji) here is given by 

A:(2,C,a;) = [ze-/VC][^e-/2C], 
and K,(z) is defined in terms of g-Pochhammer symbols [48], 

{q'z';q')oo{q^z-^;q^U 



(1.2.2) 



k{z) 



{q^z-^;q^U{q^z^;q^U' 
1 {q-'z-';q-^U{q-'z^;q-% 
L z iq-^z^;q-^U{q-^z-^;q-^)c 



\q\ < i> 
kl > 1, 



(1.2.3) 



where 



(a;6)oo = n(l-«^")- 

n=0 

This normalisation factor satisfies the functional relations 

k{z)k{1/z) = 1, 

i^{-q/z) ^ -[qz] 
k{z) [q'^/z] ' 



and has the special value 



k{1) = 1. 



(1.2.4) 



(1.2.5) 



Note that the definition (1.2.3) for k is non-analytic across the unit circle \q\ = 1. For 
these values of q there is an alternate definition of k, described in [49]. At the special 
point q = e^'^'/^, we can take k = 1, as this satisfies the functional relations (1.2.4) and 
(1.2.5). We note, however, that this choice for k is not the limit of either function in 
(1.2.3) as g ^ e^'^V^ 



Proposition 1.2.1. The Baxterised elements obey the usual Yang-Baxter and reflection 
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equations with spectral parameters: 



Ri{z)Ri+i{zw)Ri{w) = Ri+i{w)Ri{zw)Ri+i{z), 
Ko{z,ORi{zw)Ko{wX)Ri{w/z) = Ri{w/z)Ko{w,ORi{zw)Ko{z,0, (1-2.6) 
KUz, C)RL-i{zw)KLiw, QRl-i{w/z) = RL-i{w/z)kL{w, C)RL-iizw)KL{z, C). 



They furthermore satisfy the unitarity relations 

Ri{z)R,{l/z) = 1, 
ki{z,Qki{l/z,C) = l, i 



0,L. 



(1.2.7) 



Proof. These relations can be easily checked using the algebraic rules (1.1.2), or using a 
graphical notation like the one in Figure 1.6. □ 

The above Baxterised elements are special cases of i?-matrices, which can be defined 
more generally using the Hecke algebra. This is explained in more detail in Appendix A. 

We now introduce a graphical version of the Baxterised elements, using the planar 
Temperley-Lieb- Jones algebra [50] , which we will be able to use in a more general context 
than Figure 1.6. 

Definition 1.2.2. We define the i?-operator R[w, z) to be the following linear combina- 
tion of pictures: 



R{w, z) 



1 



[qz/w] 



k{w/z) \ [qw/z\ 
and graphically abbreviate R{w, z) by 

R{w, z) = 



+ 



[z/w] 
[qw/z] 



Note that we can use this picture in any orientation, as the arrows uniquely determine 
how the spectral parameters z and w enter in R. It is worthwhile to point out that 
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R{cw, cz) = R{w, z), for any constant c. We also define the boundary iiT-operators by 



fc(gw,Co,^o) 
kC^/qw, Co, ujq) 



+ 



Kl{w) 



+ 



The Baxterised elements R, Kq and Kl will be used to define the transfer matrix of 
the system. We can also write Kq as 



-qw 



qw 



(1.2.8) 



where 



Co 



+ 



which will be useful for defining the relations satisfied by Kq. 

Proposition 1.2.2. The unitarity relations (1.2.7) for R and K can he graphically de- 
picted as 



and 



(1.2.9) 






26 



In addition, the Yang-Baxter and reflection equations (1.2.6) can be written as 



and 






z w 






(1.2.10) 



We furthermore note the crossing relation satisfied by R, i.e., R{z,w) = R{—qw,z), 
which graphically reads [48] 



w = -qw 



(1.2.11) 



The above relations are straightforward to prove from the definitions, by considering 
two loops to be the same if they have the same connectivities, and replacing closed loops 
with a factor of —{q + q~^). The miitarity relation for R (1.2.9) is proved here as an 
example. 

Proof. The LHS of (1.2.9) produces four pictures, 







The first three of these have the same connectivity, so they can be grouped together with 
coefficient 



1 



-{q + Q ^)[w/z][z/w] + [qw/z][z/w] + [w/z][qz/w] 



= 0. 



k{w/z)k{z/w) \ [qz/w][qw/z] 
This leaves only the fourth picture, which is equal to the identity. □ 
Finally, we will also define slightly different versions of Kq and Kl, which will be 
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useful in proving the commutativity of the transfer matrix, 







(1.2.12) 



These i^-matriccs 



the fonowing versions of the reflection equations, 







(1.2.13) 



-yq -"^/q -yq '^/q '^q "'^Z? '^q '^/q 

These are easily proved using the crossing relation (1.2.11). 
1.2.2 Transfer matrix 

Definition 1.2.3. Using the pictorial versions of the R and K-matrices in Definition 1.2.2 
we define Sklyanin's double row transfer matrix Tl{w) = Tl{w; Zi, . . . , zl', Co, Cl) picto- 
rially in the following way (see [27,77,92]), 



Tl{w) 























































< 


































w 























This can also be written as 

Tl{w) = Try, {Koiw-^) Riw, zi) . . . R{w, zl) Kl{w) R{zl, w'^) . . . R{zi,w-^)) , 

where, in terms of pictures, the trace means that we join the two ends of the line with 
the parameter w attached to it, as shown in the above diagram. 

Proposition 1.2.3. All the possible values of w give us a commuting family of transfer 
matrices (see for example [92]), i.e., 

[Tl{v),Tl{w)]=Q, 

and hence T{w) defines an integrable lattice model. 
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Proof. First, note that with the notation (1.2.12) and the crossing relation, T(w)T(v) 
may be depicted 







































































































-v/q 















Using the unitarity relation (1.2.9) twice, and repeated application of the Yang-Baxter 
equation (1.2.10), the above becomes 























































































^ -v/q 
























-w/q 















Z7 



Now, at both boundaries, the reflection equations (1.2.13) can be applied, resulting in 























































































"-v/q 
























-w/q 















and the process involving the Yang-Baxter equation and unitarity relations can be re- 
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versed, to finally produce 











































-v/q 
























































-w/q 





























which is the graphical depiction of T{v)T{w) 



□ 



As a consequence of this commutativity, the eigenvectors of Tl do not depend on the 
spectral parameter w, but only on zi, . . . , zl, 

Tl{w;zi, ...,zl) {"^{zi, . . .,zl)) oc {"^(zi, . . .,zl)) ■ 

Following [27,29], we note that the Yang-Baxter and reflection equations (1.2.6) also 
immediately imply the following interlacing conditions of the transfer matrix with Ri, 
Kq and Kl- 

Proposition 1.2.4. 

Ri{zi/zi+i)TL{w; Zl,..., zl) = Tl{w; zi,..., Zi+i,Zi, ZL)Ri{zi/ Zi+i), 
ko{z^\Co)TL{w; Zl,..., zl) = Tl{w; z^^,Z2, Zij^oC^r', Co), (1-2-14) 
Kl{zl, Cl)Tl{w; Zl,..., zl) = Tl{w; zi,..., zl-i, z'£^)Kl{zl, Cl)- 

Pictorially, the first relation is 
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The other relations have similar depictions. 



30 



1.2.3 Hamiltonian 

Definition 1.2.4. The Hamiltonian of the TL(n) loop model with open boundaries is 
defined in terms of the logarithmic derivative of the transfer matrix with respect to w at 
the point z\ = . . . = zl = w = 1\ 



4 dw 



logr(u;;l,...,l) 
[q] IdKoiw) , dKL{w) , ^"i^dRjiw) 



w=l 

dw dw ^ dw 



+ const. 

w=l 



This can be calculated and expressed as the following operator, 

/ — sinwo \ / — sincuj^ \ 4-4 

H = co[ ^ — r - eo + CL ^ — r - cl + > ^ (-2 COS 7 - ej) + const., 

Vsm(a;o + 7) / \sm{ujL + j) J ^ 

(1.2.15) 

where cq, cl are 

4 sin 7 sin {oji + 7) 
A:(l,0,Wj) 

1.3 Specialisation n = 1 

In the next chapter we will introduce the Q'-deformed Knizhnik-Zamolodchikov equation, 
which depends on q as well as an additional parameter s. At s^q"^ = 1 the solution of this 
equation is polynomial [52,53], and in Section 2.1 we will use the interlacing conditions 
(1.2.14) to show that at = s'^ = 1 the ground state eigenvector of the TL(n) loop 
model is a solution of the qKZ equation. 

The specialisation = 1 is equivalent to setting n = 1 or 7 = 27r/3, and this 
also corresponds to anisotropy of the XXZ spin chain, A := COS7 = —1/2. In order 
to construct a special representation we will also set a;o = a;^ = —27 and 6 = 1. At 
this point we have several simplifications, the first of which is that k as given in (1.2.2) 
becomes 

kiz,C,u) = [z/qC][zC/q] =: A;(z,C). (1-3.1) 

As stated in Section 1.2.1, when q = e^'^'/^ we take the normalisation factor k{z) to 
be identically equal to 1. This choice for k causes R to be invariant under negation of its 
arguments, so the crossing relation (1.2.11) can be written without the negative sign, 

R{z,w) = R{qw,z). 

When the above specialisations are taken, the constant in the expression for the 
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Hamiltonian (1.2.15) disappears, and the rest simplifies to 

L-l 

H = co(l - eo) + cl(1 - cl) + ^(1 - e^), 

i=i 

3 

'^-i+c?+c" 

For convenience, we keep the notation and make use of the fact that q 
also see that the relations of the 2BTL algebra can now be expressed as 

ef = ei, Vi, 
ejCjiiei = Ci, 1 <i < L - 1. 

This has a one-dimensional representation p defined by 

p:ei^l, i = 0,l,...,L, (1.3.4) 

which is a quotient of the link pattern representation, as it maps every link pattern to 
the identity (this is easily seen by viewing each link pattern as a word in the e^, as 
described in Section 1.1.2). We choose the values = lol = —2^ and 6 = 1 in order 
to construct this representation. We note that p{H) = 0, hence is an eigenvalue of H 
in the link pattern representation. In fact, because the eigenvalues of Cj arc and 1, for 
co,cj^ non-negative, is the lowest eigenvalue of H and corresponds to the ground state 
of the 0(n = 1) loop model. 

Moreover, we can use a version of the Perron-Frobenius theorem to show that the 
ground state is unique. Consider the matrix H = IL — H, which has non-negative 
entries for Co,cl non- negative, and which shares eigenvectors with H. The expansion 
of H"^ includes every possible word in the generators of length m or shorter, so for 
large enough m, acting with T-L"^ on an arbitrary basis vector produces a vector with all 
entries non-zero. Thus H is irreducible, and the Perron-Frobenius theorem states that 
the eigenvector corresponding to eigenvalue L (equal to the ground state of H) is unique. 
We will be interested in the ground state eigenvector \^) as a function of the parameters 
Co and cl, 

H\^{co,cl)) = 0. 

In the same way, the one-dimensional representation indicates that the ground state 
eigenvalue of the transfer matrix is equal to 1, so the eigenvalue equation is 

Tl{w; zi,..., zl) \^l{zi, zl)) = \^l{zi, zl)) ■ (1.3.5) 

In the homogeneous limit — > 1, the transfer matrix Tl becomes the probability matrix 
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(1.3.2) 
(1.3.3) 
= 1. We 



of the stochastic raise and peel model [78,81], for which the steady state eigenvector 
(1.3.5) is unique, again by the Perron-Frobenius theorem. The eigenvalue spectrum is 
continuous, so it can be argued that there exists an open set around = . . . = = 1 for 
which the eigenvector remains unique. We will thus assume that the eigenvector remains 
unique for generic values oi zi, . . . , zl. 

The qKZ equation will allow us to obtain an explicit characterisation of | for finite 
L. We will in particular be able to derive a closed form expression for the normalisation 
Zl, which is the sum of all the components of as well as for a boundary-to-boundary 
correlation function defined in Chapter 3. In order to do so we need a recursion relation 
for I^'l), which we will discuss in Section 2.2.2. 

It is worth noting that in the limit 1 the boundary Baxterised element Kl 

maps to the identity, and this is the limit in which the two-boundary model maps to the 
one-boundary case. Similarly, in the limit eo,eL ^ 1, we obtain the trivial boundary 
case. 
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Chapter 2 



The ground state of the TL(1) 
loop model 

The aim of this chapter is to calculate the elements of the ground state eigenvector of 
the TL(1) loop model with two open boundaries, corresponding to the 2BTL algebra 
as discussed in the previous chapter. As will be shown, at = 1, the ground state 
eigenvalue equation (1.3.5) for the inhomogeneous transfer matrix of the TL(1) model 
is equivalent to a particular instance of the g-deformed Knizhnik-Zamolodchikov {qKZ) 
equation. 

The qKZ equation is a set of finite difference equations, introduced by Frenkel and 
Reshetikhin [41] in the context of the representation theory of quantum affine algebras. 
This equation depends on a number of parameters, traditionally t, q (this is not the same 
as our q), and if applicable the boundary parameters a, b, c and d. It also has two extra 
parameters, k and r, called respectively the 'rank' and 'level' of the equation. It has been 
found that when these parameters satisfy 

t''+\''-^abcd = 1, (2.0.1) 

the qKZ equation has polynomial solutions [52,53,91] (see also [38,39]). 

Interesting recent developments [53] relate polynomial solutions of the qKZ equation 
associated with Uq{slk) to the polynomial representation of the double affine Hecke alge- 
bra [51,70,88]. These solutions can be expressed in terms of Koornwinder or Macdonald 
polynomials with specialised parameters [1,39]. 

As mentioned, the ground state of the TL(1) model is a solution of the qKZ equation 
with very specialised parameters. In particular, the parameters a, b, c and d are set to 1, 
the parameter t is equal to q^ in our notation, and the parameter q is related to a new 
parameter s. We take both the rank k and the level r to be 2, and the restriction on 
the parameters (2.0.1) becomes q^s'^ = 1. Since the ground state eigenvalue equation is 
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equivalent to the qKZ equation at q"^ = 1, this condition gives us the further restriction 
s"^ = 1. As we will see, the version of the gKZ equation used here has Laurent polynomials 
as solutions. 

The ground states for TL(1) models with a variety of boundary conditions have been 
studied in the past (see for example [23,27,66, 109]), and each one is related to a version 
of the qKZ equation. In the models with zero or one open boundaries, there exists a 
highest weight vector, which is an element of the link pattern basis. The component 
of the ground state eigenvector corresponding to this basis element plays a special role 
in the solution, because it is possible to specify this component exactly, and from there 
calculate every other clement of the vector. In the case of periodic boundaries, due to 
the cyclic symmetry there is not a unique highest weight vector, but an equivalence class 
of them. Again the corresponding components of the eigenvector can be fixed, and the 
qKZ equation provides the means to then calculate the other components. 

However, in the case of two open boundaries, there is no such highest weight vector. 
This fact causes the calculation of the ground state eigenvector in the case of two open 
boundaries to be more challenging, and until recently little progress had been made. 

In this chapter we deal with this case, and generalise the results for reflecting [27] and 
mixed [109] boundary conditions. It is very hard to find a general expression for all the 
components of the two boundary eigenvector, but a certain subset of them can be found, 
and for the others certain restrictions can be made. We also calculate the normalisation 
of the eigenvector, which is the sum of all its components. This can be found explicitly 
despite the lack of an exact formula for all the elements of the eigenvector. The main 
results in this chapter were also discovered independently and through slightly different 
means by Cantini [10], in the same year. 

For small system sizes L < 2, we solve the qKZ equation explicitly for all components, 
up to an overall factor that we take to be 1. For arbitrary size, we find an explicit 
expression for two special components of the eigenvector, using recursive properties of 
the TL(1) transfer matrix. We then use the same recursive properties to find an expression 
for the overall ground state normalisation. 

With the specialisation = 1, we find that the dependence on s factors out of the 
final result. However we keep the notation for s in the expressions for the components 
with a view to generalisation. 

2.1 The g-deformed Knizhnik— Zamolodchikov equation 

As stated above, the ground state eigenvector of the TL(1) transfer matrix is a solution 
of the qKZ equation with a certain specialisation. This connection will provide the 
foundation for an explicit analysis of the ground state eigenvector for finite system size 
L. We will first describe the qKZ equation for open boundaries, and then prove the 
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equivalence with the transfer matrix eigenvalue equation (1.3.5). 

We consider a linear combination |H) of states \a) labelled by link patterns: 

\'E{zi,...,zl)) = ^^a{zi,...,ZL) \a) . 

a 

Here the sum runs over the set of link patterns of size L, and the coefficient functions 
are Laurent polynomials in the variables zi, . . . ,zl with coefficients that are functions of 
the boundary parameters and (^l as well as q and s. 

Definition 2.1.1. The g-deformed Knizhnik-Zamolodchikov equation [41,93] is a system 
of finite difference equations on the vector \E). For open boundary conditions they can 
be written as [10, 14, 27, 109],* 

Ri{zi/zi+i)\E) = 7Ti\E) , l<i<L-l, 

i^o(lM,Co)|H) =7ro|H), (2.1.1) 

Kl{sZl, sCl) \E) = TTi \E) . 

This definition uses the Baxterised elements Ri, Ko and Kl of the two-boundary 
Temperley-Lieb algebra, defined in (1.2.1). The operators Ri{zi/zi+i), Ko{l/zi,(o) and 
Ki{szl^ sQl) act on link patterns ja), whereas the operators tt^ (i = 0, . . . , L) act on the 
coefficient functions ^ only; 

TTi Zi, Zi+l, . . .) = ^(. . . , Zi+l,Zi, . . .), 1 < i < L - 1, 

TTO C{Z1,...) =^{1/Zl,...), 
TTL ^{...,Zl) =^(. ..,1/s^Zl). 

For later convenience, we note that the qKZ equations can be rewritten 

ei \E) = -Gi \E), 0<i< L, (2.1.2) 
where are the generators of the 2BTL from Definition 1.1.4, and 

/ .-.A^i/qzi+i] 1 / • / r 1 

= (TTi + 1)V^^ l<i<L-l, 

[Zi/Zi+l\ 

/ , -|^ fe(l/zi,Co) . ^ 

ao = (7ro + l) — — , (2.1.3) 

I .-.Mszl,sCl) 

a-L = -(TTL + 1) r ir 2 21 ' 

where k was defined in (1.3.1). The operators Oj {i = 0, . . . ,L) satisfy the relations of 



*We write the equations in a form used by Smirnov [93] for the type A case (see Appendix A of this 
thesis for more details on the type classifications). 
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the Hecke algebra of affine type C [23] , as well as those of the Hecke algebra of type BC. 
The Hecke algebra is described in more detail in Appendix A. 
For later convenience, we also define operators Si, such that 

Si = {q + q~^) - ai, l<i<L-l, 

sine, , ^ ^ (2.1.4) 

sm (wfe + 7) 

With our chosen specifications for q, coq, and ojl (see Section 1.3), these definitions can 
be condensed to 

Si = —1 — ai, (2.1.5) 



The Si can be explicitly written as 



\Zil Zi+l 
k{ l/sZL, sCL 

[g][s2z 



(l-TTi), l<i<L-l, 



i^y-^/ '^-^L, -^S-LJ f-, X 

^L= r„ir,2.2i (1-^l). 



2.1.1 Equivalence with the transfer matrix eigenvalue equation 

At the special values of q and s previously mentioned, the elements of the eigenvector are 
polynomials in zf, . . . ,z^, and we use this fact to construct a definition for that is 
appropriately normalised. 

Definition 2.1.2. With the specialisations listed in Section 1.3, the ground state eigen- 
vector \'^l{zi, . . . , zl)) of the transfer matrix of the TL(1) model of width L has eigen- 
value 1, and is written in the link pattern basis as 

\'^l{zi,...,zl)) = ipa{zi,...,ZL)\a) , 

aeLPL 

where the tpa are coprime polynomials. 

The basis orthogonal to the downward link patterns |q:) is given by the row vectors 
with the usual inner product (• | • ) defined by 

{j3\ a) = 5a,j3- 

Written in this basis, the left eigenvector (*| of the transfer matrix corresponding to 
eigenvalue 1 is a row vector with every element equal to 1. This is because the transfer 
matrix is a stochastic matrix, whose columns sum to 1. 
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Definition 2.1.3. The normalisation of the eigenvector is defined as 

= (*!*), 

or more exphcitly, 

Zl{zi, ...,zl) = ^ i^a{zi, ...,zl) {P\ a) 

= J2Mzi,---,zl). (2.1.6) 

a 

For the next theorem we wih also need the following lemma. 

Lemma 2.1.1. In the almost homogeneous limit zi ^ q, Z2 = ■ ■ ■ = zl = Co — Cl = ^, 
the Hamiltonian (1.3.2) can be normalised and expressed as 

H = eieo — 1. 

The ground state eigenvector of this Hamiltonian, corresponding to eigenvalue Q, is non- 
negative. 

Proof. At the point 2:2 = . . . = -zl = Co = Cl = 1) we can calculate the Hamiltonian as 
was done in Section 1.2.3, obtaining 



eieo + r^oei. 



Normalising H by ^^p^, and taking the limit as z\ — >■ g, this becomes 



H := lim H = eieo + lim \ 



21^9 [Zl\ 21^9 (1 + + Z]^ [zi] 

= eieo - 1. 

In matrix notation the non-diagonal entries of this Hamiltonian are non-negative, so 
according to the Perron-Probenius theorem the ground state eigenvector must also be 
non-negative. □ 

Theorem 2.1.2. The specialised ground state eigenvector \^) also satisfies the qKZ 
equation, Definition 2.1.1. 

The proof of this given in [22] was incomplete. We give here the full proof. 

Proof. We first act on both sides of the eigenvector equation (1.3.5) with the Baxterised 
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element R{zi/zi+i), and using the interlacing condition in (1.2.14), 



Ri{zi/zi+i) \-^L{zi,Zi+i)) = Ri{zi/ Zi+i)TL{w; Zi, Zi+i) \-^L{zi, Zi+i)) 

= Tl{w; Zi+i, Zi) Ri{zi/zi+i) Zj+i)) . 

Since the eigenvector in (1.3.5) is unique, this implies that 

Ri{,Zi/Zi+i)\^L{Zi,Zi+i)) = Pi{zi,...,ZL)\^L{Zi+uZi)) , (2.1.7) 

where Pi is some rational function. By definition, none of the elements of l^') have a 
denominator, so the denominator of ^i(. . . , Zj, 2;i+i, . . .) must be the same as the denom- 
inator of Ri{zi/zi+i), i.e., [qzi/zi+i]. We can rearrange the above, using the unitarity 
relation (1.2.7) of R, and swap the parameters Zi -H- z^+i to get 

Ri{zi/zi+i) l^^Lizi, Zi+i)) = /3i{. . . , Zi+i, Zi,..., )"^ |^'L(zi+i, Zi)) , (2.1.8) 

which leads directly to 

/3j(. . . ,Zi, Zi+i, . . ■) = Pi{. . . , Zi^i,Zi, . . . , ) ^. 

Since wc already know the denominator of Pi, the above property gives us four choices 
for the numerator, namely 

P,{z,,Zi+,) = ^^^l±lM or ±1. (2.1.9) 

[qzi/zi+i] 

In both cases the sign is fixed to +1 by setting Zi+i = zi in (2.1.8). 

Assume Pi to be the first choice in (2.1.9). In the quotient p : Cj 1 as defined in 
(1.3.4), every link pattern is projected onto the identity, and (2.1.7) becomes 

ZL{Zi,Zi+i) = ^ZL{Zi+i,Zi). 

Since is a polynomial in the z's, it must be of the form 

Zl = [qzi+i/zi]Si{zi, ...,zl), 

where Si is a polynomial symmetric in ZjjZj+i. By considering the above argument for 
all other values of i we find that Zl must be of the form 

Zl= JJ [qzj/zi]S{zi,...,ZL), 

l<i<j<L 

where 5* here is a polynomial symmetric in all its arguments. This means that the eigen- 
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vector normalisation vanishes when zi = qzj, for any j > i. However, the speciahsation 
in Lemma 2.1.1 has zi = qzj for all j > 1, and it is shown that cannot be at that 
point. Our assumption must therefore be false and we must take (3i = 1. 
Thus, we have shown that 

Ri{zi/zi+i) \^l{zi, ...,zl))= TTi \^l{zi, zl)) , l<i<L-l. 

Similarly, we can show that 

Ko{l/zi,Co) \^l{zi, zl)) = TTo \^l{zi, zl)) , 

KlIsZl, sCl) • • • , zl)) = TTL • • • , Zl)) , 

where the proof of the last equation makes use of the fact that when = 1, KL{sa, sb) = 
KL{a, h) and R{s^z, w) = R{z, w). □ 



2.2 Recursions 

We will use the qKZ equation (2.1.1) to calculate the elements of the ground state eigen- 
vector of the TL(1) loop model for size L. However this equation does not contain enough 
information to fix the elements, so we will also use a recursion that is an inherent property 
of the loop model. 

We begin by defining some maps between spaces of link patterns of different sizes. 
These maps will be very useful in constructing recursions satisfied by the transfer matrix 
and its eigenvector. 

Definition 2.2.1. For 1 < z < L + 1, let be the map that takes a link pattern of size 
L, sends site j to j + 2 for j > i, and then inserts a link from site itoi + 1, thus creating 
a link pattern of size L + 2. For example, 



f6 ■ 




I-)- 




Let ipo be the map that takes a link pattern of size L, sends site j to site j + 1 for all 
j, and inserts a ')' on the first site. 

Let ipL+i be the map that takes a link pattern of size L and inserts a '(' after the last 
site. 

We will also define similar maps on vectors in the link pattern space defined in Sec- 
tion 1.1.2. The action on the basis is a straightforward extension of the action on the 
link patterns, so we use the same notation for both maps. 
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Definition 2.2.2. Let (pi : LFl LPl^2 be the map defined by 



ipi X Ca |a) = ^Ca l^Pia) . 



a 



a 



Similarly let (po LFl ^ LP^+i and (pL+i ■ LPl LP^+i be defined by 



/3 /3 




It is worth noting that the vector resulting from the action by ipi is of length 2^+^, 
but it has only 2^ non-zero entries, which are indexed by link patterns with a small link 
from i to i + 1. Similar statements can be made for and ipL- 

2.2.1 Transfer matrix recursion 

Proposition 2.2.1. The transfer matrix satisfies the following identity for general q: 
rr I \ [q/ziw\[q^Zi/w] 

TL[w;zi,...,Zi+i =qzi,...,ZL)oipi = ^ — — - o rL_2(w;; 2:^+2, •• •)• 

[9 ZiW\[qw/ Zi\ 

The proof of this proposition is in Appendix B.l. A similar relation was proved in [27] 
and for the case of periodic boundary conditions in [29] . 

Proposition 2.2.2. Likewise, at the boundaries, and for q = e^'^^^^, the transfer matrix 



Tl{w; zi,...,ZL = Cl/q; Co, Cl) oipL = ^L° Tl-i{w; zi,..., zl-i; Co, Cl/q)- (2.2.2) 
The first of these is proved in Appendix B.2. 

2.2.2 Recursion of the eigenvector 

In order to find a recursive definition for all components of we must refer to the 

recursive property of the transfer matrix described in Proposition 2.2.1. We will suppress 
the arguments zi, . . . , zl of T and j^'i) except where detail is needed. The notation Zj 
will mean that Zj is missing from the list zi, . . . ,zl. When we specify to be a third 
root of unity, the proportionality factor in Proposition 2.2.1 becomes 1, so 



satisfies 



Tl{w; Zl =q(o,..., zl\ Co, Cl) °VQ = m° Tl-i{w; Z2,..., zl; qCo, Cl), 



(2.2.1) 



Tl{w; Zi+i = qzi) o^pi = ipiO Tl-2{w; Zi, Zi+i). 



(2.2.3) 
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Acting with both sides of this equation on the eigenvector |^x,_2(-2i, ii+i)), 
Tl{w; Zi+i = qzi) [ipi |*L-2(^i, = |*L-2(^i, ^i+i)) , 

which, by uniqueness of the eigenvector imphes 

\^L{Zi+l = qZi)) = Piizf, Zi,..., Zi, Zi+i, ...,ZL)iPi |*L-2(^i, ^i+l)) , (2.2.4) 

where pi is a proportionahty factor. This proportionahty imphes that any component 
corresponding to a Unk pattern without a small link connecting i and i + vanishes when 
Zj+i = qZi. In Section 2.4 this property of the eigenvector will be derived in another way. 
Relation (2.2.4) was already proved for subcases of the most general open boundary 
conditions in [27,109], and for periodic boundary conditions in [29]. 

Likewise, from the boundary recursions (2.2.1) and (2.2.2) of the transfer matrix we 
deduce that 

l^^iizi = gCo;Co)) = ro{z2, . . . , zl; (o) ^0 l*L-i(^i;gCo)) , 
l^^LizL = (L/q;(L)) = rL{zi, . . .,zl-i;Cl) |*l-i(^l;ClA)) , 

where vq and are proportionality factors analogous to pi. 

The above recursions for the eigenvector imply the following recursions for the com- 
ponents of the eigenvector. 

Lemma 2.2.3. 

Piizf, Zl,..., Zi, Zi+i, ...,Zl) '<Pa~'^{Zi, Zi+i), 

ro{z2,...,ZL;Co) i>0~^{zi;qCo), 

rt{zi, zl-i;Cl) V'^'H^l; CL/q), 

where a is any link pattern of length L — 2 and /5 is any link pattern of length L — 1. 

Clearly the normalisation Z given in (2.1.6) also satisfies these recursions. In Ap- 
pendix C, we show that Z is symmetric in its arguments, and use this to prove that pi is 
symmetric in all its variables except Zi, and that ro and are symmetric in all the Zi. 
It is also shown that the function pi takes the same form for each i. We henceforth drop 
the index i from p. 

*We use the term 'small' or 'little' link to mean a link connecting neighbouring sites. 
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2.3 Solutions to the qKZ equation for small L 



We now want to solve the qKZ equation to find the components of the ground state 
eigenvector. For this purpose we use the form of the equation given in (2.1.2). We will be 

looking for the lowest degree solution for which all the ipa are not identically zero. Any 
higher degree solution must be a scalar multiple of the lowest degree solution, because of 
the uniqueness of the TL(1) ground state eigenvector. 

For L = 0, the solution is trivial, but it will be needed for calculations of the recursions 
as described in Section 2.2.2. Obviously the eigenvector is of length one, i.e., a scalar. 
The definition of the eigenvector has the components being coprime, and so by analogy 
we take this scalar to be 1. 

2.3.1 Example: L = 1 

When L = 1, there are only two link patterns, denoted by ')' and '('. The ground state 
eigenvector of the TL(1) loop model in the basis of link patterns is then 

1*) = ^) \)) + ^{ |(>- 

The gj act only on link patterns, leaving the polynomials unchanged. In particular, cq acts 
on both link patterns to produce ')', and ei acts on both to produce '('. Remembering 
that the act only on polynomials, the qKZ equations can be rewritten as 

+ D) = -aoV'( 10 -ao-^) D), 

(V.( + ^))|(> = -aiV( lO-aiV^) |)>. 

Recalling that Sj = — 1 — Oj, Vi, we can match up the coefficients of each link pattern to 
get a system of four equations, 

aoV'( = 0, soV')=V'(> aiV') = 0, siip(^ = ipy 

These equations, with the definitions of Oj in (2.1.3), give us all the information we need 
to find the minimal degree solution. From the first equation, we know that 

^(1/^1, Co) V'( 

is invariant under zi — >■ 1/zi, and we can deduce from this invariance that must be of 
the form 
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where /i is a polynomial invariant under zi — >■ 1/zi. Similarly, from the third equation 
we find that ■0) must be given by 



ijj) = k{l/szi,sCL)fi{zi), 

where fi is invariant under szi — t- 1/szi. 

We find that the remaining two equations arc satisfied if fi and fi are both constants 
with fi = s^/i, and this gives us the lowest degree solution. We want the components to 
be coprime, so we choose /i to be —1,* and the solution for L = 1 is 

A = -kM)^ (2.3.1) 
= -s^k{l/szi,sCL)- 

Recursion to L = 

According to the recursions listed in Lemma 2.2.3, specialising the above solutions at 
zi = Cl/q and Zi = q(o respectively will produce the boundary proportionality factors 
ri and ro for L = 1 (recall that the solution for L = is simply 1). This gives 

rUCL) = -s^k{l/sCL,sCo), 
ro(Co) = -MCo,CL). 

2.3.2 Example: L = 2 

For L = 2, the eigenvector is of length four, and the link pattern basis consists of the 
elements! |((), |()), |)(), |))) }• As before, in the (/KZ equations the coefficients of each 
basis element can be collected together to form a system of equations. As an example, 
we take i = in (2.1.2), and obtain 

- ^ aoi/Ja |a) = X] I") 

a a 

= (V'(( + V')()|)(> + (V'() + V')))|))>- 

These equations can be written as 

= aoV'o = '^oV'((, ^2 3 3) 



* Wc choose —1 instead of 1 with the benefit of hindsight — at this stage it makes no difference to the 
solution, but this choice results in a neater expression for the final solution for general L. 
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The rest of the system for L = 2 is 



= aaV'o = «2V')), ^2 3 4) 



= aiV')) = aiV')( = ai^(( , ^ 
siV'o = '^)) + A( + ^)( • 

As in the L = 1 case, it is an easy consequence of the equation aoV'() = that if V'o / 0, it 
should contain a factor k{zi, Co)- Such vanishing conditions hold for all the components. 

Proposition 2.3.1. 

i. tpQ andtjjf^^ vanish or contain a factor k{zi, Co); ihe remainder being invariant under 
zi <H- 1/2:1. 

ii. ipQ and "0)) vanish or contain a factor k(l/sz2, sCl); the remainder being invariant 
under SZ2 -H- 1/ 52:2 . 

Hi. ■0)); and-ip^^^ vanish or contain a factor [qzi/z2], the remainder being a symmetric 
function in z\ and Z2- 

Solution 

With the known factors and symmetries from Proposition 2.3.1 above, we thus look for 
a solution of the form 



=Y[k{Zi,(o) X [qZl/z2][q/ziZ2] X f2{zi,Z2), 

1=1 

2 

1p)) = JJfe(l/'S2i,sCL) X [qZi/z2][qs'^ZiZ2] X f2{zi,Z2), 



where f2izi,Z2) is a symmetric function invariant under zi ^ 'i-/zi, and f2{zi,Z2) is 
symmetric and invariant under szi -H- 1/szi. Note that with k as defined in (1.3.1) we 
could write this as 



'^(( = n ^ Zl) X f2{zi, Z2), 

i=l 
2 

■0)) = Y[k{l/szi,sCL) X ^(1/521,52:2) X 72(21,22) 



i=l 
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The other two components may be determined from 



(2.3.6) 



We pause here to define an important polynomial, which will appear throughout the 
rest of this thesis. The symplectic character x\ of degree A = (Ai, . . . , Xl) is defined by 



deti 




deti 


L-j+l -L+j-1 





Xx{zi, ...,zl) = 

For convenience, we use the notation 

Tx{zi, . . . , zl) = x\{zl,---,zl). 

The symplectic character is symmetric, and for the degree A(^) with A 
satisfies the recursion 



(2.3.7) 



'3 



(2.3.8) 

it 



L-j 
2 



Txw{zi,...,zl) 



\zj=qzi = (-1)^ n HZi,Zk)Tx(L-2)iZi,Zj). 



(2.3.9) 



The classical character xx^ or equivalently the Schur polynomial of the type C root 
system, with degree A = A^^) appears repeatedly in related studies on loop models [27, 109] 
and symmetry classes of alternating sign matrices [71] . 

When = 1, q = e^'^'/^ we find that the solution to equations (2.3.3)-(2.3.5) can be 
given explicitly by 

f2{zi,Z2) = -r(^l,0fl){zi,Z2,CL), 
h{zi,Z2) = -T{lfi^Q){sC,0,SZl,SZ2). 

Recursion to L = 

From Condition (iii) of Proposition 2.3.1, when we set Z2 = qzi, all components vanish 
except for i/jq. From (2.3.6), 

H'^/ziXo 



k{l/szi, sCL)k{l/sz2, sCL)k{l/szi,sz2) 



X T{l,0,0)isZl, SZ2, sCo)\z2=qzi 

k{zuCofk{z,XLf, 
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where we have used the properties of k 

k{s^a, b) = s'^k{sa, sb) = k{a, 1/6) = k{l/qa, b) = k{a, b). (2.3.10) 

Since the solution for L = is simply 1, we can easily see that the proportionality factor 
p in Lemma 2.2.3 for L = 2 is 

pizi) = k{zi,Cofkizi,CLf. (2.3.11) 

Recursion to L = 1 

Similarly, when we set zi = q(o, the components ipQ and '(/'(( vanish, and V')) becomes 

V'))Ui=gCo = «^ H'^/sz2,sCL)k{Co,CL)k{Co,Z2f. 
Using our solution for tp){z2) from (2.3.1), we deduce that 

roiz2; Co) = -A;(Co, CMCo, ^2)'. (2.3.12) 
In the same way, setting 2:2 = Cl/? in V'({ gives us 

TLizi; Cl) = -s^k{l/sCL, sCo)kil/sCL, szlf. (2.3.13) 
2.3.3 Example: L = 3 

It is computationally very intensive to compute explicitly the full solution for L = 3. 
However, we can compute a solution for the subset of equations where and i{jy^^ are 
not individually determined, but only their sum is (see Appendix D for more details). 
We find 

= n ^(^J'^i) T-(1,1,0,0)(^1'^2,^;3,Cl)53(2i,^2,2;3), 

0<i<i<3 

V'))) = n H^/^^i^^^j) '^{l,l,0,0){sCo,SZi,SZ2,SZ3)gs{sZi,SZ2,SZ3), 

l<i<j<4 

with r as before, and where gs is symmetric and invariant under Zi ■<->■ 1/zi. We have 
also introduced the notation zq = Co and zl+i = Cl- Imposing the boundary recursions 
in Lemma 2.2.3 requires that the lowest degree solution is 

g3{zi,Z2,Z3) = r(^lflfi){zi,Z2,Z3). 
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Recursion to L = 2 

The above calculation also gives us 



ro{z2,Z3;Co) = ~HCo, CL)k{Co, Z2fk{Co, Z3f , 

rL{zi,Z2;CL) = -s'^k{l/ s(l, sCo)k{l/ sCl: SZlfk{l/ sCl, SZ2f . 



(2.3.14) 



Recursion to L = 1 



Computing ^'(o = •S3V'(((, and setting = qz2, we find the recursion between size L = 3 
and size L = 1: 



P{z2; Zl) = k{z2, C,Qfk{z2, C,Lfk{z2, Zi)^, 
p{zi; Zj) = k{zi, C,Qfk{zi, C,Lfk{zi, ZjY, 



(2.3.15) 



2.4 Solution for general L 

As in the case of L = 2 in Section 2.3.2, for general L we may derive factors for certain 
components. For each i from 1 to L — 1, every link pattern in the LHS of the qKZ equation 
(2.1.2) will have a small link from i to i + 1 once Cj has acted. The qKZ equation then 
says that a^'i/'a = if a does not have a small link from i to i + 1. This leads to the 
following conditions on ipa- 

Proposition 2.4.1. 

i. If a does not have a small link from the left boundary to 1, ipa vanishes or contains 
a factor k{ziXo)> the remainder being invariant under z\ -H- 1/zi. 

ii. If a does not have a small link from L to the right boundary, ipa vanishes or contains 
a factor k{l/szL, s(^l), the remainder being invariant under szl ■<->■ 1/ szl- 

Hi. If a does not have a small link from i to i + 1, ipa vanishes or contains a factor 
[qZi/zi+i], the remainder being a symmetric function in Zi and Zi+i. 

Using the above conditions, for general L the component is given by 



where //, is symmetric and invariant under Zi — > 1/zi. The majority of the factors in this 
expression are imposed by the symmetry conditions. 



(2.4.1) 



0<i<j<L 



Likewise, the component V') - ) is expressed as 



n ^(v 



SZi,SZj) fL{sZi,...,SZL), 



(2.4.2) 



l<i<j<L+l 
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where fi is symmetric and invariant under szi — t- 1/sZi. Other components may be de- 
rived from the extremal components by acting with products of Baxterised versions of the 
operators Sj, as described in Appendix A [23]. However, in the case under consideration 
it is not possible to derive every component of | ^) in this way. In Appendix D we explain 
the reasons for this in detail for the case L = 3. 

To get more information about the polynomials fi and fi, as well as about the other 
components, we use the recursive properties of the eigenvector, as we have previously for 
small L. As before, the recursions come from Lemma 2.2.3, which relates components of 
the eigenvector for L to components of the eigenvector for L — 1 and L — 2. 

2.4.1 Recursions 

We have found the proportionality factors for small system sizes, given in (2.3.2), (2.3.11), 
(2.3.12), (2.3.13), (2.3.14) and (2.3.15). We make the assumption that these factors 
continue their pattern for larger system sizes, as stated in the following conjecture. 

Conjecture 2.4.1. For general L the proportionality factors are 

p{zi;zi,...,Zi,Zi+i,...,ZL) = k{ziXofk{zi,CLf k{zi,Zj)^, 

L 

ro{z2, ...,zl; Co) = -HCo, Cl) JJ k{Co, Zif, 

1=2 

L-1 

rL{zi,...,ZL-i;CL) = -s^A;(1/sCl,<o) JJ ^(I/sCl, szj)^. 

i=l 

Using these factors, we can come up with associated recursions for the symmetric 
functions /l and /l. For example, we can calculate the size L eigenvector component 
■0^..Q using the qKZ equation and then, according to Lemma 2.2.3, setting zl = qzi-i 
will give us p{zL-i; zi, . . . , zl^2) multiplied by the size L — 2 component '4^(^.~^- Recalling 
the definition (2.1.5) of Sj, we have from the qKZ equation that 

i^f...Q = slV'(:..( , 
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and since tp^ 



(•••( 



0, it follows that 



ZL=qZL-l 



ZL=qZL-l 



n HZj,Zi)fL{zi,...,ZL) 



ZL=qZL-l 



L-2 



k{zL-l,CL)Y[k{zL-l,Zif Yi k{Zj,Zi)fL{zi,...,ZL-l,S^qZL-l) 
i=0 0<i<j<L-2 

1/ /- \2i / /- \ TJ 1 r ^2fLizi,---,ZL-l,S^qZL-l) ,L-2 

k{zL-l,Co) k{zL-l,CL) II k{ZL-l,Zi) r tp,,j . 

fJl JL-2(.Zi,...,ZL-2) ^ ^ 



Here we have used the properties of k given in (2.3.10). Prom above, the proportionality 
factor in this relation is given by p{zL-i', zi, . . . , ZL-2), so we arrive at a recursion for fi, 

L-2 

fL{zi,...,ZL-l,S^qZL-l) = k{zL-l,C,L) JJ k{zL-l, Zjf fL-2{zi, . . . , ZL-2)- 

Similarly, we can use the recursion from V'q...) to i^yT.) to find a recursion for fi,. Due 
to the symmetry properties of both these functions, the recursions can be generalised to 
arbitrary i, 

fL{zi,...,Zi,S^qZi,...,ZL) = k{Zi,CL) JJ k{Zi, Zjf fL-2{zi, . . . , Zi, Zi+i, . . . , zl), 

hizi, ...,Zi, s^qzi, ...,zl) = k{zi, Co) JJ k{zi, ZjffL-2{zi, .■■,Zi, Zi+i, zl). 



(2.4.3) 



The boundary recursions from Lemma 2.2.3 can be immediately applied to the ex- 
tremal components (2.4.1) and (2.4.2), and we find that /l and /l in addition satisfy 

L-l 

fL{zi,...,ZL-i,CL/q;Co,CL) = - n f L-i{zi, zl-V, Co, Cl/q), 

(2.4.4) 

fUqCo, Z2,..., zl; Co, Cl) = k{Co, zj) fL-i{z2, ■■■,zl; qCo, Cl), 

J=2 



where we have explicitly indicated the dependencies on Co and Cl- 
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2.4.2 Degree 

Polynomial solutions of the qWL can be labelled by their top degree /i, where is a 
partition, /xi > /U2 > . . . > /xl > 0. These solutions are of the form 

where the notation z'^ stands for the product . . . , the coefficients Cu are polynomials 
in Co ci'iid (^L, and W ■ /i denotes the orbit of under the action of the group W, defined 
as follows. The generators of W act on an L-tuple jj, by the transformations 



u 

to 

tL 



{. . . , fli, fli+i, . . .) l<i<n- 

(|Ui,/X2,...) (l-/ii,/X2,...), ifMi>0, 



where to acts as the identity if m < 0. In effect, W ■ /j. is the set of all L-tuples that can 
be obtained from /i by any combination of the above actions. 

We can use the recursions (2.4.3) and (2.4.4) to find the minimal degree of fi and fi 
for arbitrary size. The argument here is for /l, but it is easily seen that it holds for fi as 
well, and therefore that they have the same degree. In (2.4.3), consider i = 1 and denote 
the top degree of /l by i/^^^ = ■■ ■ , i'^\o, 0, . . .). The (top) degree of k{zi, Zj) is 

(1,0) in the variables zf and z?, so the degree in zf on the right hand side of (2.4.3) is 
2L — 3. Since the degree in zf on the LHS must be the same, the added degrees in zf 
and z| of /l(^i, • • • , ^l) must be greater than or equal to 2L — 3. 

In addition, by comparing degrees of any Zj in (2.4.4), it immediately follows that 
u^l^'^ is at least equal to + 1. We thus find that the following inequalities must 

hold, 



XL) 



'^2 — 



3, 



(L) . (L-l) , , 



For a minimal degree solution these inequalities become equalities, and using the solutions 
we explicitly constructed for the small system sizes L = 1, 2, 3 in Section 2.3, we find that 



L-j (j = l,...,L). 



We will write i/^^) = A^^) + A(-^+^), where A^-^) is the partition of |A(^)| = (| - l)] 
with 



L-j 



j = 1,...,L, 



(2.4.5) 
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i.e., 



A(2") = (n- l,n 1,0,0), A^^^+i) = (n, n - 1, n - 1 . . . , 1, 1, 0, 0). 

Prom the degree of k{zj,Zi) it immediately follows that the product of factors in 
the expressions for the extremal components, and V^) - )) amount to a degree of 

^(^+1) -)- A(^+^). Solutions of the qKZ equation of minimal degree, which are relevant for 
the TL(1) loop model with open boundaries, therefore have degree IJ,^^\ with 

so that 

= 2L + 1- 2j. (2.4.6) 

The total degree of these solutions is equal to \fi^^^\ = and the degree in each variable 
zf is equal to /^i = 2L — 1. 

2.4.3 Eigenvector 

By using recursion and degree properties of the general solution, we can find expressions 
for /l and /l. We emphasise again that we have taken = 1 and q = e^"^'/^. 

For L = 2 and L = 3, the solution contains a symmetric function that involves the 
symplectic character defined in (2.3.7). The solution for general L can also be expressed 
in terms of this symplectic character. It turns out that the following two functions satisfy 
the necessary recursions (2.4.3) and (2.4.4), agree with the small size solutions, and have 
the correct degree 

L(L+1) 

fL{z-i,...,ZL) = {-1) 2 T^(l+i){zi,...,Zl,(l)Tx(l){zi,...,Zl), 
fL{szi,...,szL) = i-l) 2 (s ) T^(l+i){sCo,SZi,...,SZl)t^(l){sZi,...,SZl). 

It is worthwhile noting that (2.4.3) and (2.4.4) are satisfied because of the recursion for 
the symplectic character in (2.3.9), and the specification = 1. 

2.4.4 Eigenvector normalisation 

We have derived in Appendix C, see (C.l) and (C.2), the recursions for the normalisation 
Zl as defined in (2.1.6). Using the recursion (2.3.9) for the symplectic character defined 
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in (2.3.7), we note that 

T"L+2(Co, Zi,..., ZL, Cl) Tl+i{Co, Zi, . . . , Zl) Tl+i{zi, . . . , Zl, (l) Tl{zi, ZL)\zi+i=qZi 

= k{^i, Co)^ H^i, C,Lf JJ k{zh Zj)^ tl{Co, ...,Zi, Zi+i, . . . , Cl) 

X ''"L-i(Co, ■■■,Zi, Zi+i, . . .) TL-i{. ■■ ,Zi, Zi+i, . . . , Cl) TL-2{- ■ ■ , Zi, ij+l, . . .). 

This product of four symplcctic characters therefore satisfies the same recursions as Zl 
(the boundary recursions are also easy to show). In addition, for L = 1 and L = 2 
this product is equal to Zi and Z2, respectively. Since the recursions (C.l) and (C.2) 
specify enough points to uniquely determine Zl of degree see (2.4.6), up to a factor 
independent of L, we conclude that 

Zl{zi, ...,zl) = rL+2(Co, Zl,..., ZL, Cl) tl+i(Co, zi,..., zl) 

X TL+l{zi,. . . ,zl,C,l) tl{zi, . . . ,zl). (2.4.8) 

In particular, the normalisation of the ground state of the Hamiltonian (1.3.2) is obtained 
by setting Zi = 1, and is given by 

Zl = Z2{co, cl) Zi{co) Zi{cl) Zq, (2.4.9) 

where 

Zq = tl(1, . . . , 1), 

Zl{ci) = TL+l{Ci, 1, • • • , 1), 

Z2{cq, cl) = TL+2(Co, 1, • • • , 1, Cl), 
and Cj is defined as ^ 

'^^i + C? + C^' 

as in (1.3.3). 
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Chapter 3 



The boundary-to-boundary 
correlation function 

In an attempt to rigorously prove conformal invariance of the scaling limit of critical 
two-dimensional lattice models, an interesting new class of parafermionic observables for 
Potts and loop models was recently introduced, see for example [44,86,94,96]. These 
observables are expressed in terms of the loop representation of the Fortuin-Kasteleyn 
cluster expansion of the Potts model [4, 5] , and can be shown to be discretely holomorphic 
for certain parameter values. The operators corresponding to these observables carry a 
spin conjugate to the winding angles of the loops. In this chapter we will compute an 
exact closed form expression for the expectation of one such type of observable, with spin 
one, for the two-boundary TL(1) loop model as described in Chapter 1. The remarkable 
aspect of this result is the feasibility of an exact calculation of an expectation for a finite 
size model that is not free fermionic. 




Figure 3.1: A percolating cluster with its associated hulls drawn as loops. 

We consider the TL(1) loop model as an alternative description of the two-dimensional 
critical bond percolation model [77]. In this setting the percolation model is defined on 
a finite width, vertically infinite square lattice, as in Figure 3.1, where bonds are placed 
diagonally across the faces with a certain probability. Percolation occurs when a cluster 
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of bonds stretches from the left boundary to the right. As previously mentioned, the 
loops in the TL(1) model correspond to hulls drawn around the bond clusters, so a loop 
stretching from the left to the right implies a percolating cluster. 

The expectation value calculated in this chapter can be interpreted as the density 
of percolation clusters passing in between two chosen vertices of the lattice, xi and 
X2 [12,34,66,95]. The analogous calculation for the density of clusters closed around 
a cylinder has also been performed [69]. We will sometimes refer to this density as 
the boundary-to-boundary correlation function for the two-boundary TL(1) model. Our 
calculation is valid for odd and even L, however for odd L the boundary conditions of the 
TL(1) model cause the interpretation in terms of the percolation model to break down. 

This expectation value also corresponds to the spin current in a generalised Chalker- 
Coddington network model [13, 43, 64] for the quantum spin Hall effect. The Chalker- 
Coddington network model is based on the semi-classical picture of electrons in two 
dimensions moving under the influence of a strong perpendicular magnetic field in a 
long-ranged disorder potential. The network consists of a square lattice whose edges are 
unidirectional channels and whose vertices are scattering centres. A potential correlation 
length much larger than the magnetic length leads to the formation of clusters where the 
wave function amplitude changes only slightly. These clusters correspond to the clusters 
of bonds described above. The classical limit of the generalised Chalker-Coddington 
model [64] on the square lattice is therefore described by the solvable TL(1) lattice 
model. In a further generalisation, Gruzberg et al. [43] used a pseudo spin description 
of the particle and hole states. When the SU{2) spin symmetry is observed, the full 
quantum mechanical spin current corresponds to the observable we study here. 

In this chapter we will exactly calculate the correlation function for finite system 
size L. We will first give an explicit description of this quantity, before defining it as an 
expectation value in Section 3.1.1. We then discuss the symmetries and recursion relations 
satisfied by the function, and we use these properties along with a degree argument to 
prove an explicit form for the correlation function, given in Section 3.5. Our proof is 
complete except for two conjectures, both of which we have observed to be true for small 
L, and which we expect to hold in general. 

3.1 Definition of the correlation function 

Given a configuration C of bonds on the lattice, and any two marked vertices (at positions 
xi and X2), we consider all the loops (hulls) that pass between these two vertices. We 
give a weight of to those that close on themselves or have both ends connected to the 
same boundary, and a non-zero weight to loops that connect the left boundary to the 
right. The weight given is either 1 or —1 depending on the winding of the loop, so that a 
percolating cluster with a hull that crosses back and forth between the marked vertices is 
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only counted once. The sum of the loop weights gives the number of percolating clusters 
passing in between xi and X2 in the configuration, multiplied by the sign relating to the 
winding. The average of this value over all configurations, weighted by the probability 
of obtaining each configuration, then gives the average density of percolating clusters 
between the marked points, which we refer to as the correlation function Fi^i'^a} 
This correlation function is given explicitly by 



Cer 



{xi,X2} ■ {2:1, 2:2} 

^ ^sign^ 



(3.1.1) 



Here F is the set of configurations, and the probability P{C) of a configuration C is 
the product of weights of all the individual face configurations as defined by the R and 



X-matrices in Section 1.2.1.* N, 



{xi,X2} 



C 



is the number of paths passing in between the 
points xi and X2 and running from the left to the right boundary, and signj^^'^^'^ is +1 
if xi lies in the region above the paths, and —1 if it lies below. The sign of the path in 
Figure 3.2 thus is —1. 
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Figure 3.2: Path between markers 1 and 2 (at positions x\ and X2 respectively) 
with phase —1, arising from Figure 3.1, contributing to ■ 



3.1.1 The observables 

Definition 3.1.1. We define the operator p^^'''^'^^ by marking two vertices inside the 
transfer matrix T^, from Definition 1.2.3. The locations of the markers are given by x\ 
and X2, each of which is an ordered pair (A;, y). The first coordinate, fe, can take any value 
in {1, . . . , L + 1}, corresponding to the possible positions along the width of lattice. The 
second coordinate y can take one of three values: b (bottom), m (middle), or t (top). 
The markers give weights to loops passing between them as described in the previous 
section. 

'Throughout this chapter we take q to be the third root of unity e^"^''^ and ojo = ojl = — 47r/3 as 
introduced in Section 1.3. 
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As an example of the observable we have 



^{(3,m),(3,b)} 
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The boundary-to-boundary correlation function is given by the expectation value F^^'^'^^ 
of this observable. We will now derive an expression for this expectation value. 

The probability distribution of downward link patterns at the bottom of the transfer 
matrix is given by the ground state eigenvector |*(Co) Cl! zi,. . . , zl)), calculated in the 
previous chapter, which is an element of the space LP^, as defined in Section 1.1.2. To 
describe the probabilities for the upward link patterns at the top of the transfer matrix, 
we must first define a new space, LP^, which is spanned by the upward link patterns 
of length L. The elements of this space are denoted (5|, not to be confused with the 
notation used for the left eigenvector of the transfer matrix defined in Section 2.1.1. For 
upward link patterns, we use the same shorthand notation in terms of opening and closing 
brackets as was defined in Section 1.1.2 for downward link patterns. 

Definition 3.1.2. Between LP^, and LP^ we define the LP-inner product as follows. Let 
I a) G LPi be a downward link pattern and (5| G LP^ be an upward link pattern. Then 



(a |a) 



LP 



1, V Q, a. 



This inner product can be viewed as the result obtained when an upward and a 
downward link pattern are attached to each other. All resulting closed loops and loops 
connected to the boundaries disappear, giving a weight of 1. This is best understood 
pictorially; for example, choosing a =')()())((' and a ='(())(()('5 we have 




1. 



The probabilities for the upward link patterns at the top of the transfer matrix are 
prescribed by (^'|, which is the groundstate eigenvector of the adjoint transfer matrix. 



{^TL■.= Tl{{^) = {^. 
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Lemma 3.1.1. The adjoint transfer matrix is given by 

tI{w; Co, Cl;zi,..., zl) = TL{qw''^;CL, Co; zl,.-., zi). 

Proof. This is seen by considering the result of rotating the transfer matrix by 180°. 
Each vertical pair of i?-matrices can be transformed by the crossing relation (1.2.11) into 



1/w- 



q/w ■ 



■ w/q 



The right boundary K- matrix Kl{w), when rotated by 180°, is identical to Ko{w/q) 
by (1.2.8). Similarly, Ko{l/w) rotated by 180° is identical to KL{q/w). Thus, the 
adjoint transfer matrix can be obtained by a swapping of parameters; Co Cl, -^^i -^^L, 
Z2 ■<->■ zl-1, ■ ■ ., and replacing w — >■ q/w. □ 

Lemma 3.1.2. The elements of the vector (^| are related to those of \"^) by the same 
swapping of parameters as in Lemma 3.1.1, as well as a horizontal reflection of the link 
patterns, e.g., a ='){{){)' goes to a ='{){)){' for L = 6. Thus, 

V'5(Co, Cl; Zl,..., zl) = IpaiCL, Co; ZL,..., Zl). 

The expectation value ^^^i'^^}^ ^ defined in (3.1.1), is thus given for the infinite 
lattice by 



{*|*)lp 

with defined in Definition 3.1.1. As (5| a)i,p = 1, Va, a, we have that 

(*|*)lp = 5IV^5 V'a = ^L Zl, 



(3.1.2) 



where Zl = J^a'^ct is the normalisation of the eigenvector, given by (2.4.8). Since Zl 
is symmetric in all the Zi, as well as between Co and Cl, we have Zl = Zl, so the 
denominator of (3.1.2) becomes 

additive, 

•1,2^3} 

L ~^ L —^ L 

Proof. These properties are easily seen from the definition. □ 
Because of the additivity property of Fl, it suffices to concentrate on the case where 



Lemma 3.1.3. The function Fl is antisymmetric, pjr^'^''^^^ 
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the two markers are placed on adjacent sites, so that N^^'^^^ can only be or 1. The 
markers can be separated by a horizontal lattice edge, or a vertical one as in Figure 3.2. 

There are a number of possible positions for the pair of markers in the transfer matrix. 
If the markers are vertically separated, they can be situated on the top or bottom row, 
and at any of the L + 1 positions along the width of the lattice. Recalling Definition 3.1.1, 
for the bottom row we define 

and from the expectation value yI^''\ we can calculate by 

p{ik,tUk,rn)}^^^^ . . . , zl) = yP{1/z,, . . . , 1/zl). (3.1.3) 

This fact comes directly from the definition of the transfer matrix. In addition, we will 
show (sec Section 3.2.1) that Yl is invariant under horizontal translation of the markers, 
so we sometimes drop the label k. 

Horizontally separated markers can straddle any of the L horizontal positions, and 
can be in one of three vertical positions. With the markers at the very bottom we define 

Xik) ^ ^{{k,h),{k+l,h)} 

(k) 

It will be shown in Section 3.2.2 that the expectation value Xj^ is invariant under vertical 
translation of the markers, i.e., = p{ik,™),ik+i,T[L)} _ Furthermore, 

we can write the operator X^'^ as a product, X^^ = n^x ■> where kx is the operator 
consisting of the two markers. Hence, 

L L L 

(k) 

This also implies that the expectation value Xj^ does not depend on the spectral pa- 
rameter w. 

We will at times suppress some of the arguments of polynomials. Unless other- 
wise stated, Yl = Yl{w; (q, (l', zi, . . . , zl) with a special dependence on w, and X^^ = 
zi,...,zl) with a special dependence on Zk- Where relevant, variables will 
be given explicitly, and the notation Zj in a list of arguments will mean that Zi is missing 
from the list (a convention used in the previous chapter). 

3.1.2 Relations between the observables Yl and Xl 

Relations between Xl and Yl are obtained by looking at the way the markers can move 
within configurations with fixed local orientations of paths. 
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^ ^ (k) 

Proposition 3.1.4. Yl is related to Xj^ by 

Yl{w = Zi\ zi, . . . , zl) = ^^^(-zi, • ■ • , zl), 
and X^^ is related to by 

X^l~^^\. ..,Zi = wq'^, Zi+i =w,...)oipi = ipiO Yl-2{w] ...,Zi, Zi+i, ...), 

where (pi is the map from Definition 2.2.2. 

Proof. These relations can be proved pictorially. We first prove the relation that expresses 
Yl in terms oi X^. At the specialisation w = zi, the lower face at position i only has one 
possible configuration, 



















1/W=l/Zi 
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The marker xi can be moved along the loop on this face, producing 



Zl 
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Zl 



which is just XfJ 

The relations that express Xl in terms of Yl-2 are similarly proved. For the first of 
these, the specialisation Zj+i = qzi = w imposes one possible configuration on the four 
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central faces, which acts on the small loop introduced by (fi, 



X 



(i+l) 



Zi=wq~^ 



1 


1. 
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Zi+i^qzi^w 



Zl 



Again the marker xi can be moved along the loop, and as in Proposition 2.2.1 (remem- 
bering that the chosen specialisation for q makes the proportionality factor 1), this leaves 
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Zl 



which is exactly ipi acting on Yl-2{w; ■ . ■ ,Zi, ij+i, . . .). The proof of the second relation 
is analagous. □ 



3.2 Symmetries 



(k) 

The expectation values Yl and have a number of symmetries, most arising from 

properties of the transfer matrix. In this section we will describe these symmetries. 

Lemma 3.2.1. It is easily seen that the interlacing conditions listed in (1.2.14) still hold 
to the left and right of the markers in the operators X^^ or Y^^^ , for example 

R{zi/zi+i)x'l\. . .,Zi,Zi+i, ...) = X^l'\. . .,Zi+i,Zi,.. .)R{zi/zi+i), 



for all values k ^ i,i + 1. 

Recall that the ground state eigenvector 1^*) satisfies the qKZ equation for open 
boundaries (2.1.1), explicitly 



R{zi/zi+i) \-^{. ..,Zi, Zi+i, ...)) = |*(. . . , Zi+i,Zi, . . .)) 
Ko{zT\Co)\^{zu...)) = \^iz^\...)), 
Kl{zl,Cl)\^{...,zl)) = \^{...,zI')). 



(3.2.1) 
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Similarly, the vector (\1'| satisfies 

(*(. ..,Zi, Zi+i, . . .)\R{zi+i/zi) = (*(. . . , Zi+i,Zi, ...)\, 

. . 01^0(^1, Co) = ($(^r\ • . •)!, (3-2.2) 

{^{...,ZL)\KLiz^\CL) = m---,z:^% 

These properties will be used to prove the two propositions below. 

(k) 

Proposition 3.2.2. Xj^ has the following symmetries: 

Zi+i, ...) = X^j^\. . . , Zi+i,Zi, . . .), forkj^i,i + 1, 
xP{zi,...) = xP(z^\...), forky^l, 
xi'\...,ZL)=xi'\...,z^'), fork^L. 



Therefore, x'^^^ G C[z^, . . . , i.e., regarding the variables Zk, ■ ■ ■ ,zl as complex 

(k) 

numbers, Xj^ is invariant under the action of the Weyl group of type -Bfe_i and thus is 
symmetric in ^'^Zli^j} o-s well as invariant under Zj 1/ Zj,\/1 < j < k — \ . Similarly, 
we also have that X^^^ G C[z^^j^, . . . , 2^]^^. 

Proof. For A; 7^ i, i + 1, we use the unitarity relation (1.2.7) of the i?-matrix to insert an 
identity into the numerator of X^^ ' , and by the interlacing condition in Lemma 3.2.1 and 
the qKZ equations (3.2.1) and (3.2.2), 

{^X^^^ |*)lp = {m{zi+^/zi)R{zi/zi+^)X^^^ \^)^^ 

= {^R{zi+i/zi)X^I;\zi+i,Zi)R{zi/zi+i) |*)lp 
= {^{zi+i,Zi)\x'^^\zi+i,Zi) |*(zj+i,2;i))Lp . 

(k) ■ • (k) 

Since the denominator of Xj^ ' is symmetric, we have shown that Xj^ ' is also symmetric 
under z-i -H- Zj+i. 

Similarly, for /c 7^ 1 and k ^ L respectively, 

|M/)lp = mko{zr,C^)k^{z^\Co)xf^ |*)lp 

= {^{Z^\ (zrS . . .)ko{zi\ Co) |*>LP 

= ($(zr\...,)i^f (%-\...)|^(^r\---,)>LP> 



'LP ■ 



Again, the denominator of X^^^ is invariant under z\ 1/zi and zl — > 1/-zl, so the 
result is proved. □ 
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(k) 

Proposition 3.2.3. ^^'^ following symmetries: 

yP{. ..,Zi, Zi+i, ...) = Y^\. . . , Zi+i,Zi, ...), fork^i + 1, 
Yi^\z,,...) = YP{z^\...), fork^l, 
V ..,zl) = Y^'\. . . , zl'), fork^L + l. 

Therefore we have that Y^ G C[z^, 4_ p C[z^, z^]^^ . 

Proof. The proof is exactly analagous to that of Proposition 3.2.2. □ 

(k) 

In the next section we will deduce further symmetries of Y£ \ Then in Section 3.2.2 

(k) 

these symmetries of Yl will be used to find further symmetries of ' . 
3.2.1 Symmetries of Yl 

(k) 

Proposition 3.2.4. Y^ is independent of k, and is symmetric in the combined set of 
inhomogeneities, Y^^ G Cfz^, . . . , ^;^]^^ . 

Proof. Consider an infinite lattice made up of double-row transfer matrices. Let P be 
the set of all paths p that start from the left boundary and end at the right. Since Y£ 
depends only on the spectral parameter of one double-row transfer matrix, we can assume 
that the internal spectral parameters of all transfer matrices are equal. Let us call this 
common value w. Therefore, a path is independent of its vertical position and we shall 
abuse notation to also denote by p the equivalence class of all vertical translates of a 
particular path p. The weight ujp of each path is the product of the Boltzmann weights 
associated with the local path orientations from each R and i^-matrix involved in fixing 
P- 

The set P can be split into two mutually exclusive subsets: Pt, the set of paths that 
start in the top half of the X-matrix on the left, and P\^, the set of paths that start 
in the bottom half. These sets have a 1-1 correspondence: each path p ^ P^ can be 
transformed into a path € -Pb by a vertical flip. The weights of the two paths are 
related by Up{zi, ...,zl) = • • • ' ^i})- 

Consider now Y^\ A;G{1,...,L-|-1}, and a path ^? G Pb- Because we must consider 
all of the possible vertical positions of it may have multiple contributions to Yl. The 
number of contributions rup depends on the number of times the path crosses back and 
forth over the vertical line on which the markers of Yl are placed. It can be seen that 

(k) 

when k is odd, all paths in Pb contribute positively, because each path enters Y£ from 
the left. They all contribute negatively, i.e., enter from the right, when k is even. The 
opposite is true for the paths in P^. 

For example, for A; = 1, the path pi G Pb in Figure 3.3(a) has one vertical translate 
that also winds through the two markers, hence its multiplicity m^^ = 2. The sister path 
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Pi, depicted in Figure 3.3(b), is obtained by flipping the complete configuration (but not 
the markers), and has mp-^ = 1. 





(a) path in P\, 



(b) path in Pt 



(1) 



Figure 3.3: Paths contributing to 

Since a path must cross position k from left to right exactly once more than it crosses 
from right to left, there is always one more path contributing positively than negatively. 
As a consequence of this, for p G -F^ and its sister path p e Pt, it holds that nip = rup — 1 
when k is odd and mp = nip + 1 when k is even. Consequently, Yt is given by 



Y^''\zi, . . . , zl) = (-1)'' I ^ mp up{zi,...,ZL) -'^rup Up{zi, . . . , zl) 

pePb 

As we showed in the previous section, Y^'^^ is invariant under — > l/zj, for z < A; as well 
as for i > A;, so consider the construction 

Yt\z„ ...,z,) = l (Yt\z„ ...,z,)+ Yi'\z^\ zl')) 

pen 

+ imp + {-l)'')ujp{zi, ...,ZL)-mp a;p(2;f \ . . . , z^^)^ 

= ^ (y^p{zi^,---,zl^)+'^p{zi,...,ZL)y 
p&Pb 

It thus follows that each translational equivalence class of paths p in has only one 
contribution to Yl, which is equal to the average of the weights of p and p. It is also 
clear from the above that is in fact independent of k. Since Y^ is independent of its 
horizontal position, there is no longer a restriction on the symmetries, so it is symmetric 
in Zi and invariant under Zj — )■ Vi,j. □ 



Proposition 3.2.5. As a Junction of the spectral parameter w, Yl is symmetric under 
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w q/w and a swapping of boundary parameters, 



yUw; Co, Cl; 21 . . . , zl) = YL{qw ^; Cl, Co; 2=1 • • ■ , zl)- 
Proof. Prom Lemma 3.1.1 we know that 

YLiw; Co, Cl; 21, ... , zl) = Cl, Co; zl,..., zi), 

for some A;, and by (3.1.3) this becomes 

Yl{w; Co, Cl;zi,..., zl) = YL{qw~^;CL, Co; ^/zl, lAi). 
Then use of the symmetries given in Proposition 3.2.4 completes the proof. □ 

3.2.2 Symmetries of ^ 

Proposition 3.2.6. X^^ is an element ofC[zf, . . . , Zk, ■ ■ ■ , z^]^^ , and has the property 

Proof. This comes from the symmetries of Yj, in Proposition 3.2.4, and the second rela- 
tionship in Proposition 3.1.4. □ 

Proposition 3.2.7. ' has the antisymmetry 

xi\...,z,,...) = -Xt\...,z-,\...). 

Proof. If a path p ^ Pb has a contribution of LOp{zi, . . . , zl) to x'^^\ then its sister path 
has the contribution — a;p(zf ^, . . . , z^^). As a result, X^^^ is antisymmetric when all 
Zi are sent to l/zj. Prom the symmetries already established in Proposition 3.2.6, this 
antisymmetry must be attributed to Zfc — >■ l/z^. □ 

Finally, note the additivity property around an elementary plaquette, 

X^^) _ y^'^+l) _ ^{(fc,m),(fe+l,m)} ^ ^ Q 

The invariance of Yl implies that the operator ^^('^'™)'(^+^'™)^^ with the markers in the 

^(k) 

middle of the transfer matrix, has the same expectation value as the operator Xj^ . 

3.3 Recursion relations 

The transfer matrix allows us to deduce recursion relations that the expectation values 
should satisfy, and the symmetries listed in the previous section can be used to generalise 
these relations. In this section, we will list relations that are satisfied by the expectation 
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values and Yl, in order to show that there are enough such relations, given their 

respective degrees, to determine these polynomials. Most of these relations are based 
on the recursions satisfied by the transfer matrix and the normalisation, as described in 
Section 2.2. 

Proposition 3.3.1. If the markers in the observable Fl are located sufficiently far away 
from the region where we specialise Zi, the following recursions hold, 

FL{zi+i = qzi) = FL-2izi,Zi+i), 

Fl{Co; zi = qCo) = Fi-iiqCo; h), 

Fl{Cl;zl = qQl^) = FL-iiqCj:';zL). 

For Yl, these properties always hold, as by Proposition 3.2.4 we can move the markers 
horizontally away from the specialised position. For the first recursion holds if 

k ^ i,i + 1, and the second and third hold if k ^ 1 and k ^ L respectively. 

The proof of these recursions given in [20] did not describe the recursion of the vector 
(*|. Here the full proof is given. 

Proof. Wc will prove here the first recursion. The other two are similarly proved, but 
also rely on the fact that the components of the eigenvector are invariant under (q — > 
and (l — )■ ("j^^, see (2.4.1), (2.4.2), (2.4.7), and the definition of the qKZ equation in 
Section 2.1. 

Prom Section 2.2.2, the eigenvector j^*) has a recursion at the point z^+i = qzi, 

\'^L{zi+i = qzi)) = p{zi; zi,..., Zi, Zi+i, ...,zl) ^fi \'^L-2{zi, Zi+i)) . 

From Definition 3.1.1, it is clear that the observable Fl satisfies the same recursions as 
the transfer matrix in (2.2.3), (2.2.1) and (2.2.2), 

FUzi+i = qzi) o = Q FL-2{zi, Zi+i), 

provided that the markers are not too close to i, as stated in the proposition. It remains 
to show 

($L(zi+i = qzi)\ ipi oc {^L-2izi,Zi+i)\. (3.3.1) 

We define (p\ to be a small upward link from position i to i + acting on a downward 
link pattern. We will prove here the property 

ipl l^iflizi = qZi+l)) oc \'^L-2{zi,Zi+l)) , 

and then call on Lemma 3.1.2 to prove (3.3.1). 
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The qKZ equation (2.1.2) can be written explicitly as 

a a 

and from Section 2.4, we know that the only non-zero terms on the RHS are those where 
a = ifiP for some size L — 2 link pattern /3. Thus, 

a 

where every link pattern on the RHS has a little link from i to z + 1. 

Another way of thinking of ip] is as the bottom half of an operator, which acts on 
\a) G LP/, to give \a') € LPl^2, where tpia' = e-ia. By acting on |^) with {p\ instead of 
Cj, the little link from i to i + 1 is not introduced, so in the place of the gKZ equation 
we have 

and we want to show that when specified to zi = qzi-^-i, the RHS of this equation is 
proportional to the ground state eigenvector |^'L_2(-Si, ij+i)). Recalling the definitions 
of Qi in (2.1.3), we have 

iL I ^ [Zi+l/QZi] r [qZi+i/Zi] r 

Specialising Zi = g^j+i, the second term disappears. In the first, the coefficient becomes 
1 (recall q^ = 1), and the first component recursion from Lemma 2.2.3 can be applied to 
get 

and thus, 

(p\ \^L{qZi+l,Zi+l)) = P{Zi+i; .-.^Zi, Zj+i, . . .) |^'L-2(^i, Zi+i)) . 

Using Lemma 3.1.2, this result implies (3.3.1), where the proportionality factor is p{zi). 

Finally, the proportionality factors obtained in the numerator and denominator of Fl 
cancel, resulting in the desired recursion. □ 

A description of the derivation of the bulk and left boundary recursions for the transfer 
matrix can be found in Appendix B. 
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3.3.1 Relations for Yl 

Proposition 3.3.2. Yl has the following 4L — 4 bulk recursions, 



Yl 



and for the boundaries there are 4 recursions, 



, Zj 



1 ' ' • J 1 



3 / h 



Yl{C,o; . . . , Zi, . . .) 

Yl(.Cl; ■ ■ ■ , Zi, . . . 



= YL-i{{qCo)^';...,Zi,...), 



2i=(<?Co)±i 



= YL-,i{CL/q)^';-..,Zi,...). 



Proof. By the symmetries of Yl in Proposition 3.2.4, the recursions in Proposition 3.3.1 
can be generahsed; for instance, the bulk recursions are true with Zi = qzj, not just 
Zi = qzi-i. □ 

Proposition 3.3.3. We also have the relations 

Yl{w; Co, Cl;zi,...,zl) = ±X^^\Co,Cl; zi, . . . , zl), 



Yl{w; Co, Cl;zi,. .. ,zl] 



w=qz, 



±1 



^F^l^Cl, Co; zi,...,zl). 



(k) 

Proof. These come from Proposition 3.1.4, the antisymmetry of Xj^ in Proposition 3.2.7 
and the symmetry of Yl in Proposition 3.2.5. □ 

Viewing the above relations as specialisations for Zi, and assuming that we know Xl, 
Yl-1 and Yl-2, we thus find from Proposition 3.3.2 and Proposition 3.3.3 the value of 
Yl in 4L + 4 values of Zj. 

3.3.2 Relations for 

. . (k) 

Proposition 3.3.4. Xj^ satisfies the 4L— 8 bulk recursions (we have taken i > k without 
loss of generality). 



X 



X 



(k) 



(k) 



— X^^^^i- ■ ■ , Zi, . . ■ , Zj, . . .), 



3 > k, 



Zi=q±^zf^ 



X^_2 \- ■ ■ , Zj, ■ ■ ■ , Zi, . . .), J < k. 



and 4 recursions for the boundaries, 



^L^\Co] ■ ■ ■ ,Zi,. . . 

x^^\CL;...,Zi,...) 



^= xi%{{qCof z.,...), 
^=xi%{{CL/q)^';...,Zi,...). 



^i=(?Co)='= 



Zi={CL/qy- 
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Proof. Again, the recursions in Proposition 3.3.1 can be generalised by the symmetries 

(k) 

of Xj^ in Proposition 3.2.6, provided i,j ^ k. □ 
Proposition 3.3.5. X^^ also satisfies the relations 

- YL-2{zk;Zi,Zk), 



4" 



2i = WZfc)=^^ 



provided k ^ i. 

Proof. These come from Proposition 3.1.4, and the symmetry of Yl under w ^ q/w. □ 

(k) (k) (k) 

Therefore, assuming we know Yl^2, -^L-i' know the value of Xj^ , 

viewed as a polynomial in Zi, in 4L points. 

3.4 Degrees of Yl and Xl 

We recall definition (3.1.2) of the expectation value 

Ft'^^^ = ^^{^\ Ft'^^^ I*)lp- 

Conjecture 3.4.1. The numerator of Fl contains a factor Zl, cancelling one of such 
factors in the denominator. Equivalently, F^ is of the form P/Zl, where P is a polyno- 
mial. 

For small systems we have observed this to be true, but it is still an outstanding 
problem for general L. 

(k) 

Proposition 3.4.1. Conditioned on Conjecture 3.4-1, the numerator of Xj^ must be of 
degree width AL — 2 in each variable, and the numerator of Yl must be of degree width 
4:L + 2. Thus we have enough recursion relations, listed in the previous section, to fix the 
numerators of both expectation values. 

Proof. For both expectation values X^^ and Yl, the degree width of the numerator must 
equal the degree width of the denominator. The degree width of Zl is 4L — 2 in each 

(k) 

variable Zi. By the above conjecture, this means for Xj^ ' that the numerator must also 
be of degree width 4L — 2 in each variable. 

For Yl, we must also take into account the action of the transfer matrix. The numer- 
ator and denominator of each term in the transfer matrix have degree width 4 in each 
variable. We can factor the denominator out, so the degree width of both the numerator 
and the denominator of Yl are now 4L -|- 2. □ 
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3.5 Solution 



(k) 

Before we present the exact expressions for ' and Yl, we must introduce some auxiliary 
functions. 

Definition 3.5.1. With a shght change in notation from (2.3.8), let 



tl{zi, . . . , zl) = X\iL){zl,---,zl), X 



L-j 



where the symplectic character x\ is defined in (2.3.7). 
Definition 3.5.2. We also define 

"tl+i(Co, zi,...,zl) tl+i{Cl, Zl,---, zl) 



Ul{Co, Cl; Zl,..., zl) = CL log 

where the constant cl is given by 



tl{zi, ...,zl) tl+2(Co, Cl, Zi,...,Zl) J ' 



CL = {-l)^i . 

It is interesting to note that the function ul has the form of a Toda lattice wave 
function; we refer the interested reader to [47] . 

Theorem 3.5.1. Conditioned on Conjecture 3.4-1 and Conjecture 3.5.1 (below), the 
expectation values X^^ and Yl are given by 



^L^ = Zk-^ ul{Co,Cl;zi,. . . ,zl), 
d 

Yl= w— UL+2{Co,CL;zi,...,ZL,vq~^,w) 



(3.5.1) 
(3.5.2) 



where v is a dummy variable. 

This will be proved in the next section. 

3.5.1 Proof of the main result 

It is easy to see that the degree of the numerators of (3.5.1) and (3.5.2) are 4L — 2 and 
4L + 2 respectively, as required by Proposition 3.4.1. In addition, the denominator of 
(3.5.1) is Zl, and the denominator of (3.5.2) is Zl times the denominator of the transfer 
matrix. We have also calculated by brute force the expectation values for system sizes 
L = 1,2 and 3 and verified the solution for these cases. Thus, in order to prove that 

(k) 

these expressions are the correct forms for Xj^ ' and Yl in general, we must prove that 
they satisfy all the required recursions as set out in Section 3.3. 
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Lemma 3.5.2. The function u has a list of simple recursions, based on the recursion for 
T in (2.3.9); 



Wl(Co, Cl; • ■ ■ , ZL--2, Z, qz) = UL-2{Co, Cl] ■ ■ ■ , ZL-2), 

ul{Co, Cl; zi = qCo, Z2, ■■■) = -UL-i{qCo, Cl] Z2, ■ ■ ■) - cl log [-A;(Co, Cl)] , 
ul{Co, Cl;---, zl-1, zl = q'^Cz) = -wl-i(Co, 9~^Cl; zl-i) - cl log [-k{q~'^CL, Co)] ■ 

Proof. The first of these is easily obtained by using the recursion (2.3.9) on each of the 
T- functions in the definition of , and noting that the proportionality factors on the top 
and bottom cancel exactly. The second recursion we prove now, and the third works in 
a similar way: 



ul{zi = qCo) 



Cilog 

CL log 
-CL log 



TL+l{Co, QCo, Z2,..., zl) Tl+i{Cl, QCo, Z2, - - - , Zl 
TL(gCo, Z2,..-, zl) Tl+2{Co, Cl, QCq, ^2, • ■ ■ , Zl) _ 
-1 TL-l{z2, ...,Zl) TL+l{qCo,CL,Z2, - - - , Zl 



HCo, Cl) TLiqCo, Z2,..., zl) Tl{Cl, Z2,...,Zl 

TLiqCo, Z2,-.-, zl) tl{Cl, Z2,..., Zl) 



TL-liz2, --.,Zl) TL+l{qCo, Cl, Z2, ■ - - , Zl) _ 

-UL-i{qCo, Cl; Z2, .--)-cl log [-A;(Co, Cl)] ■ 



CL log [-/c(Co,Cl)] 



□ 



(k) 

These recursions immediately imply the recursions for Yl and Xj^ ' that do not in- 
volve the marked parameter; that is, the 4L recursions in Proposition 3.3.2 for Yl-, as 
well as the 4L — 4 recursions in Proposition 3.3.4 for Xl- The remaining recursions, 
from Proposition 3.3.3 and Proposition 3.3.5, which respectively complete the degree ar- 
guments for Yl and Xj^ , are more complicated. These two recursions are addressed in 
the next section. 

Relations between Yl and ' in Proposition 3.3.3 and Proposition 3.3.5 

Conjecture 3.5.1. The proposed expression for Yl in (3.5.2) is invariant under the 
transformation w q/w. 

The proof of this in [20] was found to be flawed, and a correct proof has not yet been 
found, so we merely assume this symmetry. It can be shown explicitly that Yl satisfies 
this symmetry for small L (wc have confirmed it for L = 2 and 3), and we expect this 
symmetry to hold for arbitrary L. 

Proposition 3.5.3. The expressions (3.5.1) and (3.5.2) satisfy the first relation in 
Proposition 3.3.3, where w = zf^. Conditioned on Conjecture 3.5.1 they also satisfy 
the second relation, where w = qzf^. 

Proof. Taking expression (3.5.2) for Yl at the point w = Zi, and using the first recursion 
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in Lemma 3.5.2, 



^ ( -1 ^ 

W— UL+2{Zl,...,Zi,...,ZL,qV ,w) 



= W— UL{Zl,...,Zi,...,ZL,W) 



V=Zi 
W=Zi 



which is exactly X^f^ from (3.5.1). The other relations in Proposition 3.3.3 follow by the 
antisymmetry of (3.5.1) under Zi — >■ 1/zi and Conjecture 3.5.1. □ 



Proposition 3.5.4. The expressions (3.5.1) and (3.5.2) satisfy the relations listed in 
Proposition 3.3.5. 

Proof. We will prove the cases where z-i = Zj.q^^ and Zj = {qzk)^^, and the other two 



cases follow from symmetry under Zi 
the point Zi = ZkQ^^, 



^ We take the expression for X^^ (3.5.1) at 



Zk^ ULi...,Zi,...,Zk,...) 

dzk 

d 



Zi=Zkq 



Zk-^ UL{zi,...,Zi,...,Zk,...,ZL,Vq ^,Zk] 

dzk 



V=Zk 



which is exactly the expression for Yj, in (3.5.2) at the point w = z^. Similarly, at 
Zi = {qzk)~^, (3.5.1) becomes 



Zk^ ul{ 

OZk 



1 '^11 • • • 1 Zki 



Zi={qzk) ^ 



d 



— z 



J- ul{zi, ...,Zi,...,Zk,...,ZL,vq S ^) 



which is exactly (3.5.2) at the point w = z^ 



-1 



□ 



We expect both Conjecture 3.4.1 and Conjecture 3.5.1 to hold true for arbitrary 
system size. Under these conditions Theorem 3.5.1 is proved. 
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Chapter 4 



Separation of variables for the 
symplectic character 

In this thesis we have found exact finite size expressions for two quantities related to 
the two-boundary TL(1) loop model: the normalisation Z of the ground state eigenvec- 
tor, see (2.4.8); and the boundary-to-boundary correlation function F, see (3.5.1) and 
(3.5.2). Both of these quantities are expressed in terms of the polynomial character of 
the symplectic group, Xi defined in (2.3.7). 

For each version of the TL(1) loop model with left and right boundary conditions, 
physical quantities can be expressed in terms of this character. The function x depends on 
the bulk inhomogeneity parameters as well as, if applicable, the parameters corresponding 
to the boundaries of the lattice. Analogous calculations in the periodic version of the 
TL(1) loop model result in physical quantities that can be expressed in terms of the Schur 
function, a multivariate polynomial well-known in representation theory. 

The Schur function, corresponding to the root system of type A (in the classification 
of Dynkin diagrams) , is symmetric in all its variables, and it is the polynomial character of 
the general linear group. The symplectic character is an analogue of the Schur function 
for the root system of type C; that is, not only is it a symmetric function, but it is 
also invariant under inversion of its arguments. The boundary loop models have this 
symmetry because of the effect of the boundaries, which do not exist in the periodic case. 

In addition to appearing in solutions of statistical mechanical models, both the sym- 
plectic character and the Schur function play an important role in the representation 
theory of the classical groups [42] . They are also used as generating functions for count- 
ing problems in enumerative combinatorics; for example, each one provides a generalised 
enumeration of alternating sign matrices with various boundary conditions, as well as 
plane partitions and rhombus tilings [15,36,71]. 

For all of these applications it is important to understand the asymptotic behaviour 
of the polynomials, as the number of variables goes to infinity. In the case of the two- 
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boundary TL(1) loop model, the calculation of the asymptotic limit of the symplectic 
character will give the thermodynamic limit of the physical quantities Z and F . In 
particular we would like to analyse the effect of the two boundaries in the asymptotic limit, 
which involves setting the bulk parameters to 1 and keeping the boundary parameters 
general. 

Okounkov and Olshanski studied the asymptotic limit of type A Jack polynomials and 
type BC orthogonal polynomials with all but a finite subset of variables set to 1 [72,73]. 
For both of these problems, the authors considered polynomials with a degree growing 
linearly with the number of variables. In contrast, the symplectic character that we would 
like to analyse has a degree that is quadratic in the number of variables (see (2.4.5)) and, 
as far as we are aware, its asymptotic behaviour is an open problem. 

With the asymptotic limit described above as our motivation, in this chapter we 
study a method for the separation of variables (SoV) of symplectic characters, initiated 
by Sklyanin [92] and extended by Kuznetsov, Mangazeev and Sklyanin [56-59]. This SoV 
method makes use of the Q-operator introduced by Baxter [4] for the solution of the 8- 
vertex model. The aim of the so-called Q-operator method is to transform a chosen multi- 
variable polynomial into a product of single variable polynomials. Using this method, the 
problem of the asymptotics will be reduced to that of finding the asymptotics for each 
factorised part. In order for this approach to be useful the Q-operator method must also 
be reversible, which is one of the main technical hurdles in SoV. 

Symmetric polynomials are eigenfunctions of certain multivariate differential opera- 
tors, or Hamiltonians {i^j}, which form a quantum integrable system. The Q-operator 
for a quantum integrable system is a quantisation of the Backlund transformation for the 
corresponding classical integrable system. This connection was first found by Pasquicr 
and Gaudin [76], who discovered a correspondence between the classical Backlund trans- 
formation and the Q-operator for the periodic Toda lattice. The operator Q can be 
realised as an integral operator with a simple kernel. In the quasi-classical limit, this 
kernel turns into the generating function of the canonical transform [58, 76]. 

A number of examples are given in [59] for the application of the Q-operator method 
to systems of symmetric polynomials associated with the root system of type A. In 
particular, [59] discusses SoV for Schur functions, being a special case of Jack polynomials. 
As indicated above, this chapter aims to extend SoV using the Q-operator method to 
the irreducible character of the symplectic group, which is the Schur polynomial for 
the root system of type C . The only other study of SoV methods for root systems other 
than type A that we arc aware of is for the open quantum sl(2,R) spin chain [25]. 

The key hurdle faced when applying SoV to type A polynomials is the inversion 
of the method. The inverse of the separating operator has not been found for type 
A polynomials in general (a notable exception is the Schur function, see [59] for more 
details), and another method has to be used to reverse the process, known as the 'lifting 
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operator'. For the symplectic character, it has not been possible to find an equivalent 
lifting operator, but the inverse separating operator is easy to find, presumably because of 
the symplectic character's relationship to the Schur function. This raises some questions 
about the invertibility of the method for more general type C polynomials. 

The key ingredient of the Q-operator method is the construction of the separating 
operator S whose action on a polynomial P\ is proportional to a product of single variable 
polynomials, 

L 

(5Pa) (xi, . . . , xl) = Pa(1, • • • , 1) n '^Al^i)- (4-0.1) 

For the case of the symplectic character, as with Schur functions, S is invertible. The 
structure of the method is given in the next section. 

4.1 The Q-operator method 

We consider a quantum integrable system with a commuting set of Hamiltonians, given by 
differential operators Hi {i G {1, . . . , L}). These have eigenstates equal to the polynomial 

Py. 

HjPx{xi,...,XL) = hj{X)Px{xi,...,XL), (4.1.1) 

where the hj arc the corresponding eigenvalues, which depend on the multi-degree A = 
(Ai, . . . , Xl) of the polynomial. 

The aim of the Q-operator method is to find a version of the spectral problem (4.1.1) 
that involves a related factorised polynomial in place of Pa- Acting on both sides of 
(4.1.1) with S, 

SHjS-^ (5Pa)(xi, ...,xl) = hj{X) (5Pa) (xi, . . .,xl), 
and using (4.0.1), 

L L 
{SHjS-')l[qx{xi) = hj{X)l[qx{xi). 

i=l i=l 

The main results of this chapter are the explicit construction of the operators S and 
for the symplectic characters, as well as a 'factorised Hamiltonian' Hj, which acts in 
the same way as SHjS~^ on the factorised polynomial. Hj is not uniquely defined, and 
can be constructed using qx as shown in the next section. 
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4.1.1 Factorised Hamiltonian 

The operator Hj can be constructed from a single variable operator in the following way. 
If there exists a differential operator in Xi, denoted Wi, such that 

Wiqxixi) = 0, [Wi , qx{xk)] =0 {k i), (4.1.2) 

then we construct the operator 

Wij = Wi + hj{X). 

Noting that Wij commutes with any function of x^ where k i, it is easy to see that 

L L 

Wij n Qxi^k) = hj{X) n Qxi^k), Vi, (4.1.3) 
fc=i fc=i 

and therefore any linear combination of the form 

L L 

1=1 1=1 

will also satisfy (4.1.3). 

4.1.2 Factorisation of the separating operator S 

We will show that S can be written as a product of operators, 

5=(poQ.....Q.JL,=,^_,,=,,, (4.1.4) 

where the operator Qz^ is an integral operator that acts as 

{QziP\){xi,...,XL) = qx{zi)Px{xi,...,XL), (4.1.5) 

and po sends f{xi, . . . ,xl) — >■ /(!,...,!). Furthermore, the have the important 
properties 

Having found the operator Q^. , we introduce new operators Ai for which 

{pi-iQzi)\z.=^.=AiPi. 
Here, pj sends f{xi, . . . ,xl) f{xi, . . . ,Xj,l, . . . , 1). These new operators act as 

(A-Pa) {xi, . . . , Xi, 1, . . . , 1) = qxixi)Px{xi, Xi-i, 1 . . . , 1), (4.1.6) 
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and (4.1.4) becomes the factorisation 



S = Ai...Al. 



4.1.3 Summary 



Given a particular family of polynomials, the Q-operator method can be condensed into 
five steps: 

1. Specify Hj and hj{X) in (4.1.1), 

2. Specify Wi satisfying (4.1.2), and construct the factorised Hamiltonian, 

3. Construct S~^, 

4. Find Qz such that (4.1.5) is satisfied, 

5. Find Ai such that (4.1.6) is satisfied, and construct S. 

The inverse separating operator is given before the construction of S for two reasons; 
firstly, for the symplectic character it turns out to be simpler in form, and secondly, the 
details relate to those of the differential operators H and W used in the first two steps. 

As previously stated, this chapter focuses on the case where P\ is equal to xXi the 
symplectic character, which will be defined in Section 4.2. For this case, we closely 
follow [59], in which the factorisation of the Schur functions is treated. The technical 
detail in the symplectic case however is more elaborate. 

Section 4.2 describes the symplectic character and some of its properties. The sections 
following proceed in the order of the steps above. Section 4.5 also describes the inverse 
of the operator Sk = Ai 4fe, which satisfies 



k 



4.2 Properties of the symplectic character 

The polynomials we will consider, defined in (2.3.7), are the irreducible characters of the 
symplectic group. These form a basis of C[x^, . . . , xf^]^^, the ring of Laurent polynomials 
with complex coefficients invariant under the action of the Weyl group of type Bl, i.e., 

polynomials symmetric in Xi and invariant under Xj — t- x~^. Each character is labelled by 
a partition A = (Ai, . . . , A^), with A,j > Aj+i. We refer to them as symplectic characters 
[42,103], and for the purposes of this chapter define them as follows. 



Definition 4.2.1. 



where 



Xx{xi,...,xl) 



a^{xi, ...,xl) 



as{xi, . . . ,xl) ' 
a^{xi, ...,xl)= deti 
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6 is defined as the partition (L, . . . , L — i + 1, . . . , 2, 1), and we use fj, to denote the partition 
\ + 6, so /Uj = Xi + L — i + 1. The subscript L on the determinant denotes that the matrix 
is size L X L. 

The polynomials are Laurent polynomials in the ring C[x^, . . . that are an- 
tisymmetric in (xi, . . . ,xi^) and under Xi — t- x~^ . These antisymmetries mean that 
a^{xi, . . . ,xl) = whenever Xj = xj or Xi = 1. For the definition of Ai, we will 
need an expression for x\ when some of the arguments are set to 1, so we must find an 
alternate definition for xx in this limit. We therefore define a 'truncated' version of o^: 

Definition 4.2.2. 



a^'=)(xi,...,a;jfc_i) = det^ 



2(L-i)+l 



i<k 
i>k 

Then 



(4.2.1) 



Lemma 4.2.1. For i = k, . . . ,L, set Xi = e*^"* with ui G 

a(f)(a;i,...,a;fe_i) Ck{e)a^{xi, . . . ,xl), 
in the limit as e — >■ 0. The pref actor Ck{s) does not depend on fi. 

Proof. As an example, we first prove the above for the case where k = L. In this case 
the parameter u can be set to 1. In the Lth row of the determinant a^(a;i, . . . , e^), 
the jth element is 

e^W _ Q-^i^j = 2 sinh (e/x^) 
= 2£fij + h.o.t. 

Taking this to first order in e, we can factor 2e out of the determinant, leaving /Xj in the 
bottom row. Then 



a^ff'^xi, . . .,XL-i) = lim -^a^ixi, . . .,XL-i,e^), 
so cl{£) = l/2e. 

In the case of general k, the parameters Ui allow us to take multiple arguments to 1 
along distinct trajectories. In the rows from to L of af^{xi, . . . , e^"*^ , ■ ■ ■ , e^"^), we 
expand up to order £2(L-fc)+i^ 

+ 1)! tk (2(^ --) + !)! 
The resulting matrix can be shown to be a product of two matrices. We use 1„ for the 



1 
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identity matrix of size n, and for a matrix of zeroes whose size is determined from 
context: 



deti 



i<k 



2^m=k (2(L-m)+l)! 



i>k 



( 



det 





2i-A:+l£(-L-fe+l) 



2(£»,)^(-^-j)+l 
(2(L-j)+l)! 



,(1) 



i,j>k 



2(L-i)+l 



i<k 



i>k 



a^^\xi,...,Xk-i). 



UL=ki^iL - m) + ly. 

This leads to the result 

£-(L-fc+i)2j-|L-fc(2m + l)! , _ 
= lim ^ , ' (iV"' ^a^(xi, . . . , e^"'= , . . . , e""^), 

(Uk,-,UL) 

so the proportionality factor Cfe(e) is 

_ g-(^-^+i)^n^-i,(2m + l)! 



(4.2.2) 



Ck{£) 



□ 



We note that there is an alternative inductive proof similar to that used in [59] , which 
we give in Appendix E. 

Corollary. Because Cfc(e) is independent of ji, at the point Xk+i = . . . = = 1 we have 
an alternate definition of x, 



X\{xi, . . . ,Xfc, 1, . . . , 1) 



a}, '{xi,...,Xk) 
af^^\xi,...,Xk) 



This new definition allows us to directly calculate the symplectic character when some 
of its arguments are set to 1. 



4.2.1 Factorised forms and homogeneous identities 

This section will list some useful identities and factorised expressions for and xa- First 
recall that the denominator of a Schur function is given by the Vandermonde determinant, 
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which has the product form 

detL 



L-j 



(4.2.3) 



l<m<n<L 



The Weyl denominator formula for type C gives us an analogous identity for the sym- 
plectic character [42], 



asixi,...,XL) = Y[xi ^ Yl ( 

i=l l<i<j<L 

We can also use (4.2.3) to show 



Xi Xj J 



n {xi-xj'). (4.2.4) 



l<i<j<i 



= deiL 



2(L-j)+l 



1=1 l<i<j<L 



(4.2.5) 



This result and simple row expansion of (4.2.1) immediately leads to the following iden- 
tity, 

L 



k=l 



\ 



. j^fe l<i<j<L 



(^/*fc_^-Mfc)^ (4.2.6) 



which will be useful later. Furthermore, from (4.2.4) and (4.2.2), it is easily shown that 

L-fc k-1 



af(xi, . . . , x,_i) = n + n ^r^(^^ - 



1=0 



n ( 

l<i<j<fe-l 



i=l 



Xi X'i ) 



l<i<j<A;-l 



In particular. 



a« = n(2i + l)!, 



L-l 



which, along with (4.2.5), leads to 

xA(i,....i)=np^ n 

i=\^ \<i<j<L 

which is Weyl's dimension formula for the symplectic group [42, 103]. 



(4.2.7) 



(4.2.8) 
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4.3 Hamiltonians and eigenvalues 

We now construct the system of mutually commuting Hamiltonians Hj that satisfy the 
eigenvalue equation (4.1.1) in the symplectic case, i.e., 

HjXx{^) = hj{X)xxi^), (4.3.1) 

where we have used a bold x to refer to a list of variables xi,. . .,xl- 

Let Dx — and recall the definition of the usual elementary symmetric function 



(ij{xi, . . . , X]j) — ^ ^ Xri ■ ■ ■ Xrj ■ 

l<ri<...<rj<L 

Lemma 4.3.1. The Hamiltonians Hj, for j = 1, . . . , L, are given by 



Hjfxi^) = {as{^))-'ej{Dl,...,Dlj as{^)fx{^), 
with corresponding eigenvalues 

hjW = ej{fil,...,nl), 

recalling that fi = X + 5. 

Proof. The proof that these Hamiltonians satisfy (4.3.1) is equivalent to the proof of 

ei(-C'^i>--->^xi) a/x(x) = ej{i4,---,ljl) ^^(x). 
We set v^- = x^^ — x- and note that D^^vj = vj. Writing as 



a^ixi,...,XL)= J](-1)X'... 

aeSL 



we then have 



...... 

l<ri<...<rj<L \aeSL 



aeSL \l<ri<...<rj<L J 

Due to symmetry, the inner sum is independent of a, and equal to ej{iJ,f, . . . This 
can be taken out of the outer sum, which is equal to a^. □ 
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4.4 Differential equation for qx 

In this section we describe the single variable polynomial q\{z), and define the differential 
operator Wi that satisfies (4.1.2). 

Definition 4.4.1. The polynomial qx is given by 

/ X Xa(^, 1, • • • , 1) 
Xa(1, •••,!) 

It will be useful to introduce (j)n{zi, . . . , Zk), defined as 

0^(zi, ...,Zk) = a^J^+^\zi, Zk)/a^i}\ 

and in particular 

Z^^ — Z~^J 

where we have used the identity (4.2.6). Recalling that /j, = X + S, we can write q\{z) as 

qxiz) = M^)/Mz)- (4.4.2) 
Lemma 4.4.1. The polynomial (j^niz) satisfies the differential equation 

L 

l[{Dl-fil)Mz) = 0. (4.4.3) 

Proof. By definition, 



n=l 



so the nth term of the sum in (f)^ is reduced to by the nth factor in the product. □ 
Lemma 4.4.2. The differential equation Wq\{z) = Q is satisfied when 



n=l 

where 



=1 

(Lz^ + 2Lz -2z + L) 



- 1) 
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Proof. Equation (4.4.3) can be written as 

L 



Y[{Dl-^Jil)<t>s{z)qx{z)=G. {AAA) 



n=l 

Using (4.2.7) and (4.2.8), we obtain 

'^'^'^ = z^{2L-l)\ ' 

and applying the derivative, 

{z - l)2^-2(Lz2 + 2Lz -2z + L) 



Dz(t>5{z) 



z^{2L-l)\ 
[Lz"^ + 2Lz -2z + L) 



(z2 - 1) 

Since obeys the same product rule as the usual derivative, we have 

Mz)f{z) = Mz) [d, + ^^''^^f/i''"'^^ ) /(-) 

= Mz)Zf{z). 
Substituting this result back into (4.4.4) we obtain 

L 

0=l[{Z^-til)qx{z), 
n=l 

for zj^±l. □ 

4.5 Inverse separating operator 

The inverse of the operator S = Ai . . -Al must satisfy 

^ Xx{xi,...,xl) 



' n"^<^''= x.(i,...,i) 



for all A. 

Proposition 4.5.1. is given by the differential operator 

L 



/ 1 \ L{L — 1) -| — w- 

{S~ f) {xi,...,xl) = {-l)^^(p5{xi,...,XL)~ Kl YY4'5{xi) f{xi,...,XL), 

i=l 



83 



where 



Kl = detL 



Proof. Acting with the operator on the product of q\{xi) results, using (4.4.2), in 
Kl acting on a product of (J)f^{xi). Taking one factor of this product into each row of the 
determinant Kl, 



detL 



detL 



where we have used the exphcit expression for given in (4.4.1). This can be 

expressed as the product of two matrices, 



detL 



detL 

detL 



2(L-i)+l 



2(L-i)+l 



detL 

detL 



nm=l (/^mrin^ml/^m /^n)) 

The first determinant in the numerator is just aji\ see (4.2.5), and the second is the 
definition of c 
finally obtain 



L(L-l) Q-) „ 

definition of a^{xi, . . . ,xl)- Since the denominator is equal to (—1) 2 (a^-^) , we 



-r-r L(L-l) 

^Lll^<5(a;i)gA(a;i) = (-1) 2 0^(xi,...,xl), 



from which it immediately follows that S ^ satisfies the required property. □ 
We can also find the inverse of the operator 8^ = Ai - . -Ak, which must satisfy 



X\{xi, . . . ,Xfc, 1, . . . , 1) 
Xa(1,...,1) 



This will be useful for calculating the asymptotics of x\ with all but a few of its arguments 
set to one. 

Proposition 4.5.2. S'^^ is given by the differential operator 

k 

V) (xi, ...,xk) = i-lf-^Mxi, ■ • • , xk)-^ detfc n Mxi)fixi, Xk). 



i=l 



The power of —1 here has been corrected from that which appears in the published 
version of [21]. 
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Proof. The proof of this is equivalent to the proof of 

k 



deti 



Dlt'^\ll'^^.{x^) = (4.5.1) 



i=l 



On the LHS, we insert one factor (f)m{xi) into row i of the determinant, so that each 
element of the matrix becomes 



m 



Using the proof of Proposition 4.5.1 as a guide, we rewrite the RHS as 

fc(fc-i) 4^^a}f+^^(xi,...,a;fc) 
',(1) 

and use the product formula for the denominator while using the determinant formula 
for the numerator, to produce 



deti 



fe(fe-i) 
(-1)^ 



-^■2(L-i) + l 



i<k 



deti 



M^ 



2(L-i)+l 



(n 



k(k-l)+L(L-l) 

= (-1) ^ detL 



1^ 



2(L-i)+l 



i<k 



i>k 



detL 



2(L-j)+l 



Mfnn^i(/^f-/^n) 



k(k-l)+L(L-l) 

= (-1) ^ detL 



4L-2(i+i)+2 



(/*m-/*n) 



j>fe 



The elements in rows k to L can be written, with r] = i+ j, as 



1 



l<r<s<L 



deti 




4L-2r)" 

Mi 




" 2(L-i)- 
Ml 


j>2 . 


deti 


2(L- 


-i)" 





(4.5.2) 



which is when r] > L + 1, and 1 when rj = L + 1. This means that (4.5.2) can be 
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expressed as 



k(k-l)+L(L-l) 

(-1) 2 deti 



[*] i<k 
j<L-k 



E 



i>fe+l 
j>L-fe+l 



j>L-fe+l 



where is an entry that does not contribute to the determinant, and A is a matrix 
with I's on the backwards diagonal, O's below and '*'s above. Then the expression can 
be reduced to 



detfc 



1 rin^m 



which is exactly the LHS of (4.5.1). 



□ 



4.6 The integral operator 

In this section we will construct the operator Qz satisfying 

[QzXx] {xi,...,xl) = q\{z) x\{xi, ...,xl). (4.6.1) 

In order to construct Qz as an integral operator, we will need to define an appropriate 
domain of integration, which we will do first. The integration variables are ti, . . . ,iL-i, 
yi, . . . ,yL, and w. The integration variables interlace the x's as Xi < yi/ti < Xi+i, for 
i = l,...,L-l. 

Definition 4.6.1. For any Laurent series f{t) = ^m^j^Cmf^ with no constant term, 
i.e., Co = 0, the domain P of the integral over t is defined as 



Jp ^ Jo t r^n •''1 * 



For 1 < i < L — 1, we also define 

dt, /■*^^^+idy. 



Cm 

m 



IVi Jp H JtiXi Hi 

and we will need two domains for the integrals over yL, 

XL < < oo, and U < < xi. 

ti. .. ti-i ti . . . tL-l 

With these definitions we can now write down explicitly the operator Qz that satisfies 
(4.6.1). 
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Proposition 4.6.1. Qz is given by 



[Qzf] (xi,...,Xl) = 

1 

(t)5{z)as{xi, ...,xl) 



Jl w 



dY as{yi, . . • ,yL)/(yi, • • • 



where 



/ dY= // dyi- - // dyL_i 

rti...tL-lXi I 



Ju. 



tl...tl^_lXL 

XL TT ^itj 



dyL S iyi- wxL JJ 

\ 1=1 . 



n 



Proo/. The LHS of (4.6.1) becomes 

1 



(^5{z)a^{xx, ...,XL 
while the RHS can be written as 

a^{xi,...,XL) 
(f)s{z) as{xi, ...,xl) 

It is therefore enough to prove that 



/ / dY a^{yi,...,yL) = (t)i^{z) a,j,{xi,...,XL). (4.6.2) 

Jl ^ Jd 



To do this, we first use the following three ingredients: 
i. We expand a^(y) along the bottom row, producing 

L 



^(-l)^+'(y£'-y^^'-)detz.-i 



r=l 



(4.6.3) 



For each term in this sum, the integrals over can be performed easily, resulting 
in 



-nf) -(in-,) 



(4.6.4) 



ii. The integral over w becomes elementary: 



dw 
Jl w 



{w'^'" + w"'^'-) = — {z^'- - z' 

fir 



(4.6.5) 
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iii. Combining (4.6.3) and (4.6.4) we find the further simphfication 



fir 



-FT — 

XL l\ Xitf 











detL-1 





L 

- JJx^-'^'- detL_i 



after which the remaining integrals over Yi can be inserted into each row of the 
determinants, e.g., 



dY^-- dyL_idetL_i 



^2Atr^yMj-Mr _ y-fJ,j-fJ,r\ 



= detL_i 



:det 



L-l 



JJ' j - IJ-r 
(-1)^-1 



-fjj-flr 



'flj—flr\ 



detL- 



-flj. 



Using these ingredients, we finally find the following expression for the LHS of (4.6.2), 



E 



i-iy-' fn^r det^-1 l^^i^Zi - - ^i'^^i^i' - 



Kl=l 
L 



n^r'^detz.-! x^^^-x,^D-:^{^ -x,^^) 



(4.6.6) 



At this point it has become clear that if (— 1)*'^^ times the difference of products is 
independent of r (so that we can factor it out of the sum), the remaining sum will be 
precisely equal to 4>fj,(z). Hence, it remains to show that for each r, 



a^(x) = i-iy-' detL-i k^^f C 



n ^7'^ deti_i a;f;, (xfj, - x^^ ) - (xf - x;^^ ) 
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This can be achieved by increasing the size of the matrices in the determinants by 
one row and one column, at i = 1 and j = r, cancelhng the factor of (— 1)*^'^^. The 
new column has entries of except at i = 1, which is 1, and the new row has entries 
x^'^''{xi^ — x^^^) for the first determinant and x'^''{xi^ — x^'^^) for the second. We then 
use row reduction, adding the 1st row to the 2nd row, then the 2nd to the 3rd, and so 
on. For ease of display, we place the additional column on the end, producing an extra 
factor of (—1)^"'' that will be cancelled at the very end when we permute the column 
back to j = r. Thus we obtain 



(-l)^"MlI^^detL 



"''{Xi' -x^ 







1 









1=1 



where [1] denotes a column of I's. The factors outside the determinants are then inserted 
into their respective rows, resulting in 



(-1) 



L-r 















- deti 























-fir 



These two determinants can be combined by simply performing the subtraction in the 
last column. Permuting this column back to j = r, this finally results in 



deti 



a^,{xi,...,XL). 



□ 



4.7 The integral operator A]^ 

The next and final step is to obtain the operator Ak, which satisfies 

{pk-iQzfiy^))\^^^^ = Akpkfi^), 

where 

Pkf{xi,...,XL) = f{xi,...,Xk,l,...,l)- 

Proposition 4.7.1. Relation (4.7.1) is satisfied by 



Ak = lim 



<t>s{xk)a\ '{xi,...,Xk. 



(4.7.1) 
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where 



/ dY(^) = // dYi--- dYk-i 



rti-.tk-ixi / ge^^2^^~ 



+ dyk S iyk - we^Yi 

Jti ...tt._i \ 7_i 



tfxi 



and 



/ti...t/c-i e« y yi ^ 

dtk-i dyk-i 



J Jvk-i Jp ^k-l Jtk-iXk-i Vk-l 

We remind the reader that the domain of integration P is given in Definition 4-6.1. 

Proof. The x\ form a linear basis of C[y^, . . . ,y^]^^, so it suffices to show that (4.7.1) 
holds for f = ^^d the general case will follow by linearity. We thus want to prove 
that AkPkXxi^) = qx{xk)pk-iX\{^): i-e-, 

lim ^ r d y(^) a(f+i)(j/i, ...,yk) = M^khPix^, . . .,Xk-i), (4.7.2) 

using the fact that af^'^^^ pkXx = ci^ix'^^''- 

This proof is very similar to that of Proposition 4.6.1, but with a few added subtleties. 
The differences in the proof are outlined here, and we give the details in Appendix F. 

As in step i. of the proof of Proposition 4.6.1, the determinant a^i^~^^\xi, . . . ,Xk) is 
expanded along a row, but this time along the kih row instead of the last. Also, the 
integrals over y^ take the place of that over in this step. The simplification in step 
iii. is also very similar. The products outside the determinant now run from 1 to A: — 1, 
and contain an extra factor of c'^^'*''. The k — 1 integrals over Yi, . . . , Y^-i are inserted 
into the first k — 1 rows of each determinant. 

After these integrals are performed, the sizes of the matrices are increased as before. 
This time we choose the entries of the new row proportional to sinh(e//j) so that the 
matrices simplify after row reduction. After this, the prefactors in front of the two 
determinants can be cancelled, and the determinants can be combined. 

The final step is to approximate to first order in e, and then perform row reduction 

on the rows from k to L, in order to remove the dependence on fj,r from each row. The 

(k) 

result is a factor (xi, . . . , Xk-i), which can be factored out of the remaining sum over 
r. This last sum is simply equal to (t>n{xk), and the proof is complete. □ 
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Conclusion 



Two physical quantities of interest for the Temperley-Lieb loop model at the special point 
n = 1 have been calculated in this thesis, the groundstate eigenvector and its normalisa- 
tion in Chapter 2 and the boundary-to-boundary correlation function in Chapter 3. In 
addition, a key step has been laid down for the asymptotic analysis of these results in 
Chapter 4. 

We have shown in Chapter 2 that the ground state eigenvector equation for the 
TL(1) loop model with open boundaries is equivalent to the g-deformed Knizhnik- 
Zamolodchikov equation. We have used this connection to obtain an explicit description 
of key components of the ground state |^) for finite system sizes and without resorting 
to the Bcthc ansatz. The derivation of the normalisation Z of the TL(1) ground state 
is presented in Section 2.4.4:, with the result expressed as a product of four symplectic 
characters x (also known as Schur polynomials of the type C root system). Our results 
for two open boundaries contain as special cases those of reflecting and mixed boundary 
conditions, which have been considered before [23,27,109]. The general result for the 
open boundary case has proved to be less tractable than those corresponding to these 
other boundary conditions, as it is not possible to express the eigenvector components 
as factorised operators acting on a 'highest weight' component. Though we were able 
to find symmetries and relations that the eigenvector components must satisfy, in some 
cases the relations did not contain enough information to specify the components exactly. 

Our result for the normalisation leads in Chapter 3 to an explicit expression for a 
finite size correlation function F of the open TL(1) loop model. In a two-dimensional 
percolation model on a finite width lattice, this quantity is the average number of per- 
colating clusters passing in between two chosen points. This observable also falls within 
a class of parafermionic correlation functions that can be shown to be discretely holo- 
morphic, and hence are precursors to analytic correlators in the conformally invariant 
continuum limit. We stress furthermore that our result is exact for systems of finite 
width, and not asymptotic, which is unusual for systems that are not free fermionic. Our 
main results (3.5.1) and (3.5.2) are expressed in terms of the same symplectic character 
X as appeared in the normalisation. Analogous results have been obtained for cylindrical 
boundary conditions for both site and bond percolation [69]. 
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Naturally we desire to study the asymptotic limit of these results, as this limit corre- 
sponds to the thermodynamic limit of the system. In the case of the correlation function, 
we wish to use the asymptotic analysis to check our results against the continuum limit, 
which was calculated non-rigorously using conformal field theory by Cardy [12]. To calcu- 
late the asymptotic limit of Z and F, an asymptotic analysis of the symplectic character 
X must first be conducted. We are particularly interested in the role the two boundaries 
play in the limit as L goes to infinity, so we would especially like to study the asymptotics 
of the symplectic character when all zi are set to 1 but Co and Cl are kept general. 

In Chapter 4 we have set out the Q-operator method of separation of variables for 
the symplectic character, as initiated by Sklyanin [92] and as applied to type A Jack 
polynomials in [56,59]. The steps of the separation method arc laid out in Section 4.1. 
The separation of variables for the symplectic character follows a process similar to that 
of the type A Schur polynomial (more commonly known as the Schur function). However, 
some of the intermediate steps are technically much more involved. In particular, the Q- 
and .4-operators contain a double integration in each variable, whereas the corresponding 
operators for the Schur function contain only one. This is related to the fact that the 
Hamiltonian for the symplectic character contains double derivatives. 

The problem of asymptotics for the symplectic character can now be regarded as 
a problem of asymptotics for the separated polynomial, which is a product of single- 
variable polynomials q\. As mentioned earlier, we are particularly interested in the 
asymptotic limit of x\ when all but k variables are set to 1, and the operator of 
Proposition 4.5.2 and its inverse are useful here. The differential equation satisfied by 

given in Section 4.4, is expected to be helpful in determining the asymptotics of q\. 
This is certainly worthy of further investigation. 
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Outlook 



The XXZ quantum spin chain with open boundaries at both ends is closely related to the 
TL(1) loop model at the anisotropy point A = —1/2, as discussed in the Introduction. 
We hope and expect that the results in this thesis will lead to the calculation of finite 
size correlation functions for the XXZ spin chain. The symplectic character can also be 
interpreted as a Toda lattice wave function, an observation which hints at an interest- 
ing link between the quantum integrable TL(1) model and classical integrable models. 
Other incarnations of the TL(1) model to which our results may be applied include the 
conformally invariant stochastic raise and peel model [78,81] and supersymmetric lattice 
models with boundaries [19,107]. 

In the case of the stochastic raise and peel model, the ground state is a stationary 
state distribution. In this setting an expression for the normalisation is important, so 
that the stationary state is a properly normalised probability distribution. 

One way forward would be to express certain linear combinations of the components 
of the ground state eigenvector in terms of multiple contour integrals such as was done for 
reflecting boundary conditions [31,110]. However, the success of this approach seems to be 
related to the existence of factorised expressions for the ground state components, which 
exist for reflecting and mixed boundary conditions [23]. In the case of open boundaries, 
the lack of a pseudovacuum leads to a lack of such convenient factorised expressions, 
so we anticipate some fundamental difficulties with this approach. We further hope to 
make a connection between our solutions and Macdonald-Koornwinder polynomials of 
type (C^, C) for specialised parameters [51], as well as with those in the form of Jackson 
integrals for qKZ equations on tensor product spaces, see [65,99] for the case of type A. 

There is another, independent reason for computing the normalisation, which is in 
the context of the Razumov-Stroganov conjecture [3, 17,84]. The original version of this 
conjecture was proved in [11], but the more general versions remain unproved. This 
conjecture states that there is an intriguing relation between the TL(1) ground state 
I^'l) and the combinatorics of a fully packed loop (FPL) model on flnite geometries, 
as well as other combinatorial objects such as alternating sign matrices and symmetric 
plane partitions [9]. In particular it states that the ground state normalisation is equal 
to the statistical mechanical partition function of an FPL model on a finite patch of 
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the square lattice. Such a hnk has been made for the ground state of the model with 
identified open boundaries, see [24], but for the model considered here, which has two 
genuinely open boundaries, it is not known which FPL geometry (if any) gives rise to a 
partition function equal to the normalisation of as given in (2.4.9). Understanding 
the underlying combinatorics for the general case of two boundaries will therefore lead to 
a deepening of our understanding of the Razumov-Stroganov correspondence, as well as 
to possible generalisations of symmetric plane partitions and alternating sign matrices. 

The method of separation of variables for the symplectic character is a small part of a 
much larger theory of polynomials. An obvious extension of this work is to generalise the 
method to Jack polynomials of type C and BC, or even Koornwinder polynomials. It is 
not expected that the inverse separating operator will be easy to construct in the more 
general case, as it is still an unsolved problem in the case of type A Jack polynomials. 
However, another method for reversing the SoV process was used in [56], using a so-called 
lifting operator, and it is hoped that an analogue of this will prove useful for the more 
general C and/or BC polynomials. 
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Appendix A 

The Hecke algebra 



/2-matrices are fundamental objects in the study of integrable lattice models [2]. In a 
lattice model of width L, Ri describes the possible states at site i. The i?-matrix defined 
in (1.2.1) is just one example of a more general i?-matrix, 

Ri{z) = a{z) + b{z)Ti, 1 < i < L - 1, 

where the Tj are algebraic objects, and a(z) and b{z) are coefficient functions. The 
Tj satisfy certain relations such that when a{z) and b{z) are chosen appropriately, the 
Yang-Baxter equation, 

Ri{z)Ri+i{zw)Ri{w) = Ri+i{w)Ri{zw)Ri+i{z), 

is satisfied. 

If the model in question is periodic, the site i = L is identified with i = 1, and the 
Yang-Baxter equation is also valid for Rl-i and ^i. 

If the model has non-diagonal boundary conditions, iC-matrices are introduced to 
describe the possible states at one or both of the boundaries. A more general form of a 
-fC-matrix is 

Ki{z) = c{z) + d{z)Ti, i = 0,L. 

Again, c{z) and d(z) are coefficient functions and the algebraic object Tj is chosen so 
that the reflection equation, 

Ko{z)Ri{zw)Ko{w)Ri{w/z) = Ri{w/z)Ko{w)Ri{zw)Ko{z), for i = 0, 
Kl{z)Rl-i{zw)Kl{w)Rl~i{w / z) = Rl-i{w / z)Kl{w)Rl-i{zw)Kl{z) , for i = L, 

is satisfied. 

The relations satisfied by the Tj are those of the Heckc algebra [63,89]. The Hecke 
algebra is related to the symmetric group, and it has a number of versions relating to the 
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Dynkin diagram 'type' classifications. 

Definition A.l. The Hecke algebra of type is generated by {Tj \ l<i<L — 1}, 
and defined by the relations 

{Ti - q){Ti + g"^) = 0, TiTi+iTi = Tj+iTiTj+i, 

with the Ti commuting otherwise. The latter of these relations is commonly known as 
the braid relation. 

The most relevant presentation for the Temperley-Lieb loop model is given by the 
simple transformation 

ei=Ti-q, (A.l) 

so that 

ef = -{q + q~^)ei, ejCj+iei - = ej+ieiej+i - Cj+i. 

By setting eiCj+iej — = 0, we obtain the 5/2 quotient known as the (zero-boundary) 
Temperley-Lieb algebra given in Definition 1.1.1, with q = e^'^. A representation of the 
Temperley-Lieb algebra on the tensor product space (gj^^C^ is 



and is used to build the i?-matrix of the six vertex model, Rq^. 

By taking slk quotients of the Hecke algebra, one can build modules of rank k — 1. 
Such higher rank modules arc discussed in [32,91]. For example, one can construct 
a representation on the tensor product space ^^^iC^, which can be used to build an 
i2-matrix generalising Rq^ [6,54]. 

In addition, systems with boundaries can be dealt with by including additional gen- 
erators in the algebra, with boundary-type braid relations. We give one example here. 

Definition A. 2. The Hecke algebra of type Bl has the generators and relations of 
Definition A.l, along with the additional generator Tq that satisfies 

(To - e''^°){To - e- ^'^o) = 0, TqTiToTi = TiToTiTq. 

The presentation defined by (A.l) can be extended to the type B case by defining 

2 i sin {ujo + 7) ' 





1 q 

g-^ 1 





1 



(g)L-i-l 
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so that 

2 — sin cjo 

Co = ; — rco, eoeieoei - cqCi = eicoeieo - eieo- 

sm (wo + 7) 

Again, the s/2 quotient is obtained by setting both sides of this last relation to 0, and 
this produces the one-boundary Tempcrlcy-Lieb algebra from Definition 1.1.3. 

Similarly, the Hecke algebra of type BCl is defined by the generators and relations 
of type Bl, along with the additional generator T^, which satisfies similar relations to 
To, with Tl-1 taking the place of Ti. However, the type BCl algebra is no longer finite- 
dimensional, as elements of the form {TqTi . . . T^)^ cannot be reduced. The presentation 
defined in (A.l) can again be extended, and the SI2 quotient produces the two-boundary 
Temperley-Lieb algebra from Definition 1.1.4 [18]. 

Finally, it is worthwhile to point out that the operators a, and Sj defined in (2.1.3) 
and (2.1.4) are generators of the Hecke algebra of type EC [23]. They are related to the 
Ti by 

ai = q-Ti, Si = q-^+Ti, l<i<L- 1, 

flfc = ^. ■ . — ^, Sk = . — r, k = 0,L, 

2 1 sm [LOk + 7) 2 1 sm [ujk + 7) 

and satisfy 

ttiSi = SiQi = 0, < z < L, 

along with the relations (we list them here for the Oj; the relations satisfied by the Si are 
identical) 

= + 'r^)o-ii 1 < ^ < - 1, 



2 sm Wfe 

sm iiOk+l) 



2 OllitV/c 

Ofc = ZTTF-. — TTT^fe' k = l},L, 



and 



aiai+iUi — tti — aj+iajaj+i — Cj+i, 1 < i < L — \, 

a^aiaoai — a^ai = oiaooiao ~ Q-iao, 
QLaL-iaLaL-i - aiai-i = ai-iaiaL-iaL - ai-iai- 

For the zero-boundary and one-boundary q¥JL equations (corresponding to the type 
A and type B Hecke algebras respectively), the highest weight element of the solution 
vector can be found explicitly by fixing its degree and considering its symmetries. All 
the remaining elements can then be written in terms of operators acting on the highest 
weight element. These operators are Baxterised versions of the Hecke algebra generators 
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Si [23], and have the form 



Because the type BC Hecke algebra is infinite-dimensional and lacks a highest weight 
element, this approach only partially works for the qKZ equation with two boundaries. 
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Appendix B 

Proofs of the transfer matrix 
recursions 



B.l Proof of Proposition 2.2.1: The bulk recursion 

In the following we will use the shorthand notation 

"^"^ M' '^'^ [qzY 

We first define a'^ = </?taL-2 to be the link pattern of length L with a small link connecting 
sites i and i + 1 inserted into the link pattern aL-2- Restricting our focus to the action 
of the transfer matrix on the sites i and i + 1, we have 



TL{w;zi+i = qzi) [u'l) 



















w 


-4 












W 






















— ► 
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As each ii-operator consists of two terms, the action of on sites i and i + 1 produces 
sixteen terms, 



a{ziw)a{qziw)a{w / Zi)a{w / qzi) 



+ 



a{ziw)a{qziw)a{w/zi)b{w/qzi) ^^^^ 



Some of these pictures are equivalent with respect to their external connectivities. In total 
there are five different kinds of connectivities. For instance, one of the connectivities has 



a{ziw)b{qziw)a{w / Zi)b{w / qzi) 



a{ziw)b{qziw)b{w / Zi)b{w/ qzi) 



a{ziw)b{qziw)a{w / Zi)a{w / qzi) 





4 






■ w ■ 






























; 

r 













+ 



+ 



The closed loop in the first diagram is erased at the expense of a factor — (g + g"^), after 
which the coefficients of the three diagrams sum to 0. Using the fact that a{qu)a{u) + 
b{qu)b{u) — {q + q~^)a{qu)b{u) = 0, it is easy to show that this happens for four of the 
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five total kinds of connectivities, and that we are then left with 
TL{w;zi+i = qzi) [a'^) 

























a{ziw)b{qziw)b{w / Zi)a{w / qzi) * 
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= a{ziw)b{qziw)b{w/zi)a{w/qzi) ipi 
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— >■ 



[q/ziw] [q^Zi/w] 



cpi Tl-2{w\ Zi, Zi+l) \aL-2) 



[q'^Ziw] [qw/zi 

B.2 The left boundary recursion (2.2.1) 

Here we will describe the recursion of the transfer matrix at the left boundary. The 
recursion at the right boundary is similar. 

Setting zi = q(o, the entries of the eigenvector are zero unless they correspond to 
link patterns with a little link from position 1 to the left boundary. We have the transfer 
matrix acting on this little link, 



21=9C0 



There are 8 possible configurations in the above picture, which can be grouped into 
3 kinds of connectivities. The following is a summary of the connectivities along with 
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their weights: 




+ 



qw/qCo]J \k{'w/q, Co) Kw/q.Co) 



[q^Cqw] \[qw/qC,Q 
^ [1/gCo^'] [q^Co/w] k{q/w,Co) 
[q^Cow] [qw/qCo]k{w/q,Co) 

[1/Cow] , [ l/qCow] [q'^Co/w] [qwCo/q] [qw/ gCo] 
[q^Cow] [w/Co] [q^/wCo] [q^Co/ w] 

[wCo] 



+ 



[q'^CoM 

[1/CoH ^ 
[q'^Cow] ' [q^Cow] 



+ 



0, 



[l/qCow] [qCo/w] f k{q/wXo) ^ [l/g][g'>' 
[q^Qow] [qw/qQa] \k{w/q,Co) k{w/q,Co) 

^ [qC.olw][i/q(,Qw] 

^ k{q/w,qCo) 
k{w/q,qCo)' 



i^/qCow] [q\o/w] [l/q][qyw^] 
[q'^Cow] [qw/qCo] k{w/q,Co) 

_ [l/qCowmo/w][l/q][qyw^] 

[q'^Cow] [w/ Co] [q'^Co/w] [q'^/wCo] 

[w/CoWCow] 
[l/q][qVw^] 



k{w/q, gCo) 



The second and third of these can be seen as the two contributions to a left boundary 
E'-matrix, with an inserted loop connecting to the left boundary. The weights agree, and 
the whole sum can therefore be written as ipQ o Kq{w, qQo), in a system of size L — 1. This 
implies 



= (pooTL-i{qC,o;z2,. . .)■ 
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Appendix C 



Symmetries of Z and the 
proportionality factors p^, rg, and 

Recall the ground state eigenvector normalisation, 

= (*!*), 

which is equal to the sum over the components of \^). Acting with the left eigenvector 
of the transfer matrix, (^'|, on both sides of the gKZ equation (2.1.1), and using the fact 
that TTj commutes with (^'|, 

a 

Since ei\d) = \a') for some link pattern a', and (*| a) = 1 for all a, this becomes 

a 

Similar arguments can be made to show that ttqZl = Zl and ttlZ^ = Z^. We therefore 
know that Zl remains unchanged under any permutation of the variables Zi. Recalling 
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that 



we have 



ZL{zi+i = qzi) =Pi{zi; 



Zi, Zi+i, . . .)ZL-2{Zi, Zi+l)- 



(C.l) 



Since Z is symmetric in the Zj, pi becomes pj under the transformation Zi -H- Zj. We 
henceforth drop the index i from pj. Further, since the LHS is symmetric in everything 
except Zi, as is p{zi; . . .) must also be symmetric in all its variables except for Zj. 

Completely analogously we can derive boundary recursions for the normalisation, 
which are 

Zl{zi = qCo,..., zl; Co, Cl) = ro{z2, ..■,zl; Co)Zl-i{z2, ■■■,zl; gCo, Cl), 

(C2) 

Zl{zi,...,zl = ClA;Co,Cl) = rL{zi, . . . , ZL-i; Cl)Zl-i{zi, . . . , ZL-i; Co, CL/q)- 

As before, this shows that the proportionality factors ro and ri must be symmetric in 
all their variables except for (q and (l respectively. 
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Appendix D 

Solving the gKZ equation for L = 3 



Here we will solve the qWL equation (2.1.2) with L = 3. For this example we will take 
q = e^^^i/^ and = 1, but we will leave b generic. We use the notation 

V'2 = i'io V'6 = '0)0 

V'4 = V'()) V8 = V'))), 

and we recall the definition of Si in (2.1.5) 

Si = -1- ai. 

Considering in turn each i and a, the qWL equation, 

XI V'a (e? |a)) = - X (diV'a) la) , 

a a 

implies the following system of 32 equations: 

) a = l,...,4 (D.l) 



SlV'4 = V'6 + '08 



1, 2, 5, . . . , < 



(D.2) 
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a2V'a = a = 1, 3, 4, 5, 7, 8 

S2lp6 = 1p4 + 1p5 + + 'ipS (D-3) 

a3^2a = q; = 1,...,4 

S3V'2a-l =1p2a a = 1, 3 (D.4) 

S3tp2a-l = bi)2a a = 2,4: 

The relations where the action of the projector Oj is zero entails certain symmetry re- 
strictions on the components. For instance, 

3 

l<j<j<3 i=l 

3 

l<i<j<3 i=l 

where the fa{zi, Z2, z^) are symmetric functions invariant under Zi — t- 1/zi. 

Using the system of equations we may obtain the components ipa in terms of V'l and 
■08- Assuming wc know V-'i, we find ■i/'2 = •ssV'i , then tp^ = S21P2 — i^i, and tp^ = b~^ssips. 
Then we can find tps from 

■08 = 821^6 - V'4 - V'5 - (D-5) 
by using iJjq = siip4 — tps and applying S3 on both sides of (D.5), 

-1p8 = S3S2{sitp4 - tps) + V'4 - ■06 - btps, 

which implies 

(6 - l)lp8 = {S3S2S1 - Si + 1)'04- 

The expressions for 1^2, tp3, '04 as well as ipg can be neatly rewritten in a factorised form 
as in [23]. However, if 6 = 1, the component ips cannot be determined this way. 

In a similar way, given ■08 we can find tp^, then ipQ, and ip2- Since we can express 
■04 and ip2 in two different ways, these components have to satisfy certain consistency 
conditions. Now we can find the remaining two components, 

(6 - 1)'05 = Sl1p3 -1pl-1p2- S2lp6 + ■04 + ^08, 
■07 = •Sl'03 -tpl-1p2- btp5- 

Again, if 6 = 1, -05 and ipj cannot be found separately in this way. However, their sum 
can be determined. Assuming we find an expression for one of these components (say, by 



115 



solving so^5 = ipi for ^5) that satisfies the appropriate degree, the boundary recursions 
in Lemma 2.2.3, and the symmetries imposed by (D.2) and (D.3), the entire system can 
be shown to be consistent. 

At 6 = 1, as L becomes larger there are more and more components like -05 and 
-07, which cannot be calculated directly from the extremal components. However, as 
above, certain sums of these components can be found, and consistency conditions can 
be checked. Using the recursions discussed in Section 2.2.2, recursive conditions can also 
be put on these components. 
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Appendix E 

Inductive proof of Lemma 4.2.1 



Proof. This proof shows that o[f^(xi, . . . , x^-i) is proportional to a^J^'^^\xi, . . . , x^) in 
the limit as Xfe — > 1, and relies on the proof for A; = L given in Section 4.2. 

We take lim£_^o ojf^^^ {xi, ■ ■ ■ , Xfe-i, e^). In the fcth row of the determinant (given in 
(4.2.1)), the jth element is 



e^.. - e-^.- = 2 sinh {s„j ) = 2 X; • 



\2n+l 



We use row reduction with the rows below the fcth to remove the terms up to and including 
n = L — (k + l). The remainder of the series is 

^2(L-^ + l)l +''-°-*- 

SO al!'\xi, Xk-i) = lini£->o ^2£2(^-feHi ' "m""^^^ (^i ' ■ • • ' ^k-i,e^)- Then by induction we 
have that 

a(f)(xi,...,a;fc_i) oc lim a^(xi, . . . , ojl). □ 
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Appendix F 

Proof of Proposition 4.7.1 



Proof. The proof of equation (4.7.2) is similar to that of Proposition 4.6.1: 

i. The first step is to expand the determinant a/j over the kth row, giving 



E(-i)'+^(yr-yr'")det^-i 



flj -flj 
Vi -Vi 



2(L-i)+l 



i<k 



i>k 



(F.l) 



We can perform the integrals over in (4.7.2), and, as before, the dependence on 
w factors out. 

ii. The integral over w has the same evaluation as (4.6.5). 

iii. The remaining factor is combined with the determinant. 



k-l 



- e 



fc-i 

n 



yi 



1=1 'I 



tfxi 



deti-i 



i<k 



2(L-i)-l 



i<k 



k<i<L 



2(L-i)-l 



k<i<L 



k-l 



- i<k 



2(L-i)-l 



k<i<L 



The k — 1 integrals over in (4.7.2) can be inserted into the first k — 1 rows of the 
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determinants, and evaluated as before. The two determinants then become 



k-1 

n 

1=1 



deti-i 



i<A;-l 



^(2eT^'^'- sinh(6/x,-) - 4-1 - xZi)) 



-IX 



2(L-i)-l 



fe-l<i<L 



X deti-i 



2eT-M. sinh - x^^^ (x^i, - x.^O 



i<fc-l 



k-l<i<L 



3+r 



Again, as in (4.6.6), we will extract the ingredients needed for </>^(xjfc) and show 
that the rest is independent of r so that it can be factored out of the sum in (F.l). 

As in Section 4.6, wc increase the size of each matrix by adding a row at i = L and a 
column at j = r, introducing a factor of (—1)^"'"''. The new column has entries equal to 
except at i = L, which equals 1, and the new row has entries of 2e^^'^^ sinh(£/Xj). We 
depict column r at the end: 

k-l 



1=1 



X det^ 



i+l -^i+l 



2 e^-M- sinh {s^ij) - xl'^ (x^i^ - x^^{ 



i<k—l,j^r 



k—l<i<L,j^r 



[2e=F^Mrsinh (enj)] 








1 



Now we use row reduction, subtracting the Lth row from the A;th and then adding 
the (fc — l)th to the A;th, the (fc — 2)th to the (A; — l)th, etc. We then multiply the first 
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k rows by —x^^^' , and the Lth by e^'^'^'' 
(-l)^-MetL ' 



i<k,jy^r 



[2 sinh {eiij)\ 



k<i<L,j^r 



h 





i<k 



At this point we are able to combine the two determinants by subtracting the rth column 
of the second from the rth column of the first, obtaining 



deti 



i<k,j^r 



{x: 



i<k 



[2sinh(e//j)] 



k<i<L,jjtr 







2 sinh (e/ir) 



Now we approximate to first order in e, resulting in 2e/Xj in the bottom row. We can 
factor out 2e, and then use row reduction once again: To each row in turn from i = L — 1 
to i = k we add /x^ times the row below, and we are finally left with 



2£ deti 



i<k 



2{L-i)+i 



i>k 



2ea{f)(xi,...,a;fc_i; 



This then factors out of the sum in (F.l), which is equal to (j)^{xk). Putting everything 
together, the factors of (—1) cancel out, as does the factor of 2e, and we are left with 

a'j^^xi, . . .,Xk-i)4>f,{xk), 



which is the RHS of (4.7.2). 



□ 
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